Capítulo 1 


LA COMPLEJIDAD DE LOS ALGORITMOS 


1.1 INTRODUCCIÓN 


En un sentido amplio, dado un problema y un dispositivo donde resolverlo, es 
necesario proporcionar un método preciso que lo resuelva, adecuado al dispositivo. 
A tal método lo denominamos algoritmo. 


En el presente texto nos vamos a centrar en dos aspectos muy importantes de los 
algoritmos, como son su diseño y el estudio de su eficiencia. 


El primero se refiere a la búsqueda de métodos o procedimientos, secuencias 
finitas de instrucciones adecuadas al dispositivo que disponemos, que permitan 
resolver el problema. Por otra parte, el segundo nos permite medir de alguna forma 
el coste (en tiempo y recursos) que consume un algoritmo para encontrar la 
solución y nos ofrece la posibilidad de comparar distintos algoritmos que resuelven 
un mismo problema. 


Este capítulo está dedicado al segundo de estos aspectos: la eficiencia. En 
cuanto a las técnicas de diseño, que corresponden a los patrones fundamentales 
sobre los que se construyen los algoritmos que resuelven un gran número de 
problemas, se estudiarán en los siguientes capítulos. 


1.2 EFICIENCIA Y COMPLEJIDAD 


Una vez dispongamos de un algoritmo que funciona correctamente, es necesario 
definir criterios para medir su rendimiento o comportamiento. Estos criterios se 
centran principalmente en su simplicidad y en el uso eficiente de los recursos. 


A menudo se piensa que un algoritmo sencillo no es muy eficiente. Sin 
embargo, la sencillez es una característica muy interesante a la hora de diseñar un 
algoritmo, pues facilita su verificación, el estudio de su eficiencia y su 
mantenimiento. De ahí que muchas veces prime la simplicidad y legibilidad del 
código frente a alternativas más crípticas y eficientes del algoritmo. Este hecho se 
pondrá de manifiesto en varios de los ejemplos mostrados a lo largo de este libro, 
en donde profundizaremos más en este compromiso. 


Respecto al uso eficiente de los recursos, éste suele medirse en función de dos 
parámetros: el espacio, es decir, memoria que utiliza, y el tiempo, lo que tarda en 
ejecutarse. Ambos representan los costes que supone encontrar la solución al 
problema planteado mediante un algoritmo. Dichos parámetros van a servir además 
para comparar algoritmos entre sí, permitiendo determinar el más adecuado de 
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entre varios que solucionan un mismo problema. En este capítulo nos centraremos 
solamente en la eficiencia temporal. 


El tiempo de ejecución de un algoritmo va a depender de diversos factores 
como son: los datos de entrada que le suministremos, la calidad del código 
generado por el compilador para crear el programa objeto, la naturaleza y rapidez 
de las instrucciones máquina del procesador concreto que ejecute el programa, y la 
complejidad intrínseca del algoritmo. Hay dos estudios posibles sobre el tiempo: 


1. Uno que proporciona una medida teórica (a priori), que consiste en obtener una 
función que acote (por arriba o por abajo) el tiempo de ejecución del algoritmo 
para unos valores de entrada dados. 


2. Y otro que ofrece una medida real (a posteriori), consistente en medir el tiempo 
de ejecución del algoritmo para unos valores de entrada dados y en un 
ordenador concreto. 


Ambas medidas son importantes puesto que, si bien la primera nos ofrece 
estimaciones del comportamiento de los algoritmos de forma independiente del 
ordenador en donde serán implementados y sin necesidad de ejecutarlos, la 
segunda representa las medidas reales del comportamiento del algoritmo. Estas 
medidas son funciones temporales de los datos de entrada. 


Entendemos por tamaño de la entrada el número de componentes sobre los que 
se va a ejecutar el algoritmo. Por ejemplo, la dimensión del vector a ordenar o el 
tamaño de las matrices a multiplicar. 


La unidad de tiempo a la que debe hacer referencia estas medidas de eficiencia 
no puede ser expresada en segundos o en otra unidad de tiempo concreta, pues no 
existe un ordenador estándar al que puedan hacer referencia todas las medidas. 
Denotaremos por 7(n) el tiempo de ejecución de un algoritmo para una entrada de 
tamaño n. 


Teóricamente T(n) debe indicar el número de instrucciones ejecutadas por un 
ordenador idealizado. Debemos buscar por tanto medidas simples y abstractas, 
independientes del ordenador a utilizar. Para ello es necesario acotar de alguna 
forma la diferencia que se puede producir entre distintas implementaciones de un 
mismo algoritmo, ya sea del mismo código ejecutado por dos máquinas de distinta 
velocidad, como de dos códigos que implementen el mismo método. Esta 
diferencia es la que acota el siguiente principio: 


Principio de Invarianza 


Dado un algoritmo y dos implementaciones suyas /, e >, que tardan T,¡(n) y Tan) 
segundos respectivamente, el Principio de Invarianza afirma que existe una 
constante real c > 0 y un número natural ny tales que para todo n > ny se verifica 
que T¡(n) <€ cTXn). 


Es decir, el tiempo de ejecución de dos implementaciones distintas de un 
algoritmo dado no va a diferir más que en una constante multiplicativa. 
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Con esto podemos definir sin problemas que un algoritmo tarda un tiempo del 
orden de T(n) si existen una constante real c > 0 y una implementación / del 
algoritmo que tarda menos que c7(n), para todo n tamaño de la entrada. 


Dos factores a tener muy en cuenta son la constante multiplicativa y el ny para 
los que se verifican las condiciones, pues si bien a priori un algoritmo de orden 
cuadrático es mejor que uno de orden cúbico, en el caso de tener dos algoritmos 
cuyos tiempos de ejecución son 10%? y Sn* el primero sólo será mejor que el 
segundo para tamaños de la entrada superiores a 200.000. 


También es importante hacer notar que el comportamiento de un algoritmo 
puede cambiar notablemente para diferentes entradas (por ejemplo, lo ordenados 
que se encuentren ya los datos a ordenar). De hecho, para muchos programas el 
tiempo de ejecución es en realidad una función de la entrada específica, y no sólo 
del tamaño de ésta. Así suelen estudiarse tres casos para un mismo algoritmo: caso 
peor, caso mejor y caso medio. 


El caso mejor corresponde a la traza (secuencia de sentencias) del algoritmo que 
realiza menos instrucciones. Análogamente, el caso peor corresponde a la traza del 
algoritmo que realiza más instrucciones. Respecto al caso medio, corresponde a la 
traza del algoritmo que realiza un número de instrucciones igual a la esperanza 
matemática de la variable aleatoria definida por todas las posibles trazas del 
algoritmo para un tamaño de la entrada dado, con las probabilidades de que éstas 
ocurran para esa entrada. 


Es muy importante destacar que esos casos corresponden a un tamaño de la 
entrada dado, puesto que es un error común confundir el caso mejor con el que 
menos instrucciones realiza en cualquier caso, y por lo tanto contabilizar las 
instrucciones que hace para n= 1. 


A la hora de medir el tiempo, siempre lo haremos en función del número de 
operaciones elementales que realiza dicho algoritmo, entendiendo por operaciones 
elementales (en adelante OE) aquellas que el ordenador realiza en tiempo acotado 
por una constante. Así, consideraremos OE las operaciones aritméticas básicas, 
asignaciones a variables de tipo predefinido por el compilador, los saltos (llamadas 
a funciones y procedimientos, retorno desde ellos, etc.), las comparaciones lógicas 
y el acceso a estructuras indexadas básicas, como son los vectores y matrices. Cada 
una de ellas contabilizará como 1 OE. 


Resumiendo, el tiempo de ejecución de un algoritmo va a ser una función que 
mide el número de operaciones elementales que realiza el algoritmo para un 
tamaño de entrada dado. 


En general, es posible realizar el estudio de la complejidad de un algoritmo sólo 
en base a un conjunto reducido de sentencias, aquellas que caracterizan que el 
algoritmo sea lento o rápido en el sentido que nos interesa. También es posible 
distinguir entre los tiempos de ejecución de las diferentes operaciones elementales, 
lo cual es necesario a veces por las características específicas del ordenador (por 
ejemplo, se podría considerar que las operaciones + y + presentan complejidades 
diferentes debido a su implementación). Sin embargo, en este texto tendremos en 
cuenta, a menos que se indique lo contrario, todas las operaciones elementales del 
lenguaje, y supondremos que sus tiempos de ejecución son todos iguales. 


Para hacer un estudio del tiempo de ejecución de un algoritmo para los tres 
casos citados comenzaremos con un ejemplo concreto. Supongamos entonces que 
disponemos de la definición de los siguientes tipos y constantes: 
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CONST n =...; (* num. maximo de elementos de un vector *); 
TYPE vector = ARRAY [1..n] OF INTEGER; 


y de un algoritmo cuya implementación en Modula-2 es: 


PROCEDURE Buscar(VAR a:vector;c: INTEGER) : CARDINAL; 
VAR j:CARDINAL; 


BEGIN 
j:=1; (* 1 +) 
WHILE (a[j]<c) AND (¿<n) DO (* 2 *) 
j:=3+1 (* 3 x) 
END; (* 4 x) 
IF a[j]=c THEN (x* 5 x) 
RETURN j (* 6 *) 
ELSE RETURN O (* 7 x) 
END (* 8 x) 
END Buscar; 


Para determinar el tiempo de ejecución, calcularemos primero el número de 


operaciones elementales (OE) que se realizan: 


En la línea (1) se ejecuta 1 OE (una asignación). 


En la línea (2) se efectúa la condición del bucle, con un total de 4 OE (dos 
comparaciones, un acceso al vector, y un AND). 


La línea (3) está compuesta por un incremento y una asignación (2 OE). 
La línea (5) está formada por una condición y un acceso al vector (2 OB). 
La línea (6) contiene un RETURN (1 OE) si la condición se cumple. 


La línea (7) contiene un RETURN (1 OE), cuando la condición del /F anterior 
es falsa. 


Obsérvese cómo no se contabiliza la copia del vector a la pila de ejecución del 


programa, pues se pasa por referencia y no por valor (está declarado como un 
argumento VAR, aunque no se modifique dentro de la función). En caso de pasarlo 
por valor, necesitaríamos tener en cuenta el coste que esto supone (un incremento 
de n OE). Con esto: 


En el caso mejor para el algoritmo, se efectuará la línea (1) y de la línea (2) sólo 
la primera mitad de la condición, que supone 2 OE (suponemos que las 
expresiones se evalúan de izquierda a derecha, y con “cortocircuito”, es decir, 
una expresión lógica deja de ser evaluada en el momento que se conoce su 
valor, aunque no hayan sido evaluados todos sus términos). Tras ellas la función 
acaba ejecutando las líneas (5) a (7). En consecuencia, 7(n)=1+2+3=6, 


En el caso peor, se efectúa la línea (1), el bucle se repite n—-1 veces hasta que se 
cumple la segunda condición, después se efectúa la condición de la línea (5) y la 
función acaba al ejecutarse la línea (7). Cada iteración del bucle está compuesta 
por las líneas (2) y (3), junto con una ejecución adicional de la línea (2) que es 
la que ocasiona la salida del bucle. Por tanto 
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T(n) = (0) +4) +2+1=6n+2. 


¡=1 


e En el caso medio, el bucle se ejecutará un número de veces entre O y n—l, y 
vamos a suponer que cada una de ellas tiene la misma probabilidad de suceder. 
Como existen n posibilidades (puede que el número buscado no esté) 
suponemos a priori que son equiprobables y por tanto cada una tendrá una 
probabilidad asociada de 1/n. Con esto, el número medio de veces que se 
efectuará el bucle es de 


Tenemos pues que 
(m1) /2 Y 
T(n)=1+ (CD í4+12))+42) 1241=30+3. 
i=l 


Es importante observar que no es necesario conocer el propósito del algoritmo 
para analizar su tiempo de ejecución y determinar sus casos mejor, peor y medio, 
sino que basta con estudiar su código. Suele ser un error muy frecuente el 
determinar tales casos basándose sólo en la funcionalidad para la que el algoritmo 
fue concebido, olvidando que es el código implementado el que los determina. 


En este caso, un examen más detallado de la función (¡y no de su nombre!) nos 
muestra que tras su ejecución, la función devuelve la posición de un entero dado c 
dentro de un vector ordenado de enteros, devolviendo O si el elemento no está en el 
vector. Lo que acabamos de probar es que su caso mejor se da cuando el elemento 
está en la primera posición del vector. El caso peor se produce cuando el elemento 
no está en el vector, y el caso medio ocurre cuando consideramos equiprobables 
cada una de las posiciones en las que puede encontrarse el elemento dentro del 
vector (incluyendo la posición especial 0, que indica que el elemento a buscar no se 
encuentra en el vector). 


1.2.1 Reglas generales para el cálculo del número de OE 


La siguiente lista presenta un conjunto de reglas generales para el cálculo del 
número de OE, siempre considerando el peor caso. Estas reglas definen el número 
de OE de cada estructura básica del lenguaje, por lo que el número de OE de un 
algoritmo puede hacerse por inducción sobre ellas. 


e Vamos a considerar que el tiempo de una OE es, por definición, de orden 1. La 
constante c que menciona el Principio de Invarianza dependerá de la 
implementación particular, pero nosotros supondremos que vale 1. 


e El tiempo de ejecución de una secuencia consecutiva de instrucciones se calcula 
sumando los tiempos de ejecución de cada una de las instrucciones. 
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e El tiempo de ejecución de la sentencia “CASE C OF v1:S1lv2:S2]...|lvn:Sn 
END;” es 7 = T(C) + max (7(S¡),7(82),...,1(S,)). Obsérvese que 7(C) incluye el 
tiempo de comparación con v,, va ,..., Vp. 


e El tiempo de ejecución de la sentencia “IF C THEN Si ELSE $52 END;” es 
T= T(C) + max (1(S1),1(S2)). 

e El tiempo de ejecución de un bucle de sentencias “WHILE C DO S END;” es 
T= T(C) + (n” iteraciones)*(7(S) + 7(C)). Obsérvese que tanto 7(C) como 7(5) 
pueden variar en cada iteración, y por tanto habrá que tenerlo en cuenta para su 
cálculo. 


e Para calcular el tiempo de ejecución del resto de sentencias iterativas (FOR, 
REPEAT, LOOP) basta expresarlas como un bucle WHILE. A modo de ejemplo, 
el tiempo de ejecución del bucle: 


FOR i:=1 TO n DO 
S 
END; 


puede ser calculado a partir del bucle equivalente: 


i:=1; 

WHILE i<=n DO 
S; INC(i) 

END; 


e El tiempo de ejecución de una llamada a un procedimiento o función 
F(P,, P»..., P,) es 1 (por la llamada), más el tiempo de evaluación de los 
parámetros Py, P,,..., P,, más el tiempo que tarde en ejecutarse FF, esto es, 
T=1+7(P,) + T(P,) + ... + T(P,) + T(£). No contabilizamos la copia de los 
argumentos a la pila de ejecución, salvo que se trate de estructuras complejas 
(registros o vectores) que se pasan por valor. En este caso contabilizaremos 
tantas OE como valores simples contenga la estructura. El paso de parámetros 
por referencia, por tratarse simplemente de punteros, no contabiliza tampoco. 


e El tiempo de ejecución de las llamadas a procedimientos recursivos va a dar 
lugar a ecuaciones en recurrencia, que veremos posteriormente. 


e También es necesario tener en cuenta, cuando el compilador las incorpore, las 
optimizaciones del código y la forma de evaluación de las expresiones, que 
pueden ocasionar “cortocircuitos” o realizarse de forma “perezosa” (lazy). En el 
presente trabajo supondremos que no se realizan optimizaciones, que existe el 
cortocircuito y que no existe evaluación perezosa. 


1.3 COTAS DE COMPLEJIDAD. MEDIDAS ASINTÓTICAS 


Una vez vista la forma de calcular el tiempo de ejecución 7' de un algoritmo, 
nuestro propósito es intentar clasificar dichas funciones de forma que podamos 
compararlas. Para ello, vamos a definir clases de equivalencia, correspondientes a 
las funciones que “crecen de la misma forma”. 
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En las siguientes definiciones N denotará el conjunto de los números naturales 
y R el de los reales. 
1.3.1 Cota Superior. Notación O 


Dada una función f, queremos estudiar aquellas funciones g que a lo sumo crecen 
tan deprisa como f. Al conjunto de tales funciones se le llama cota superior de f y lo 
denominamos O(f). Conociendo la cota superior de un algoritmo podemos asegurar 
que, en ningún caso, el tiempo empleado será de un orden superior al de la cota. 


Definición 1.1 


Sea f: N>[0,00). Se define el conjunto de funciones de orden O (Omicron) de f 
como: 


O(f)= (e: N>[0,00) | JceR, c>0, ImeN + g(n) <c Ain) Va 2 m0). 
Diremos que una función £: N >[0,00) es de orden O de f'si £ e O(f). 
Intuitivamente, t e O(f) indica que £ está acotada superiormente por algún 


múltiplo de f. Normalmente estaremos interesados en la menor función f' tal que £ 
pertenezca a O(f). 


En el ejemplo del algoritmo Buscar analizado anteriormente obtenemos que su 
tiempo de ejecución en el mejor caso es O(1), mientras que sus tiempos de 
ejecución para los casos peor y medio son O(n). 


Propiedades de O 


Veamos las propiedades de la cota superior. La demostración de todas ellas se 
obtiene aplicando la definición 1.1. 


1. Para cualquier función f'se tiene que fe O(f). 

. FeO0(g) > O(f) < O(g). 

. O(f) =0(g) S fe O(2) y g €0(/. 

. SifeO(g) y € e€0(h) > f € O(h). 

. SifeO(g) y fe 0(M) > f € O(min(g,h)). 

. Regla de la suma: Si f¡ €O(g) y £ € 0(4) > f + f € O(max(g,h)). 

. Regla del producto: Si f¡ eO(g) y f4 € 0(M) > f1:f € O(g:h). 

8. Si existe lim 4) 
">= g(n 


a) Sik 0 y k < oo entonces O(f) = O(g). 


b) Si k = 0 entonces f eO(g), es decir, O(f) C O(g), pero sin embargo se 
verifica que g ¿ O(f). 


A DD Uu Bb UN 


= k, dependiendo de los valores que tome k obtenemos: 
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Obsérvese la importancia que tiene el que exista tal límite, pues si no existiese 
(o fuera infinito) no podría realizarse tal afirmación, como veremos en la 
resolución de los problemas de este capítulo. 


De las propiedades anteriores se deduce que la relación —o, definida por f—o g sl 
y sólo si O(f) = O(g), es una relación de equivalencia. Siempre escogeremos el 
representante más sencillo para cada clase; así los órdenes de complejidad 
constante serán expresados por O(1), los lineales por O(n), etc. 


1.3.2 Cota Inferior. Notación Q 


Dada una función f, queremos estudiar aquellas funciones g que a lo sumo crecen 
tan lentamente como f. Al conjunto de tales funciones se le llama cota inferior de f 
y lo denominamos Q(f). Conociendo la cota inferior de un algoritmo podemos 
asegurar que, en ningún caso, el tiempo empleado será de un orden inferior al de la 
cota. 


Definición 1.2 


Sea f: N>[0,00). Se define el conjunto de funciones de orden Q (Omega) de f 
como: 


Q(A = fe:N >[0,0) | 3ceR , c>0, Ime Ne g(n) >c fam) Vn >n0). 


Diremos que una función £:N >[0,00) es de orden Q de f'si £ e Q(f). 


Intuitivamente, £ e Q(f) indica que £ está acotada inferiormente por algún 
múltiplo de f. Normalmente estaremos interesados en la mayor función f tal que £ 
pertenezca a Q(f), a la que denominaremos su cota inferior. 


Obtener buenas cotas inferiores es en general muy difícil, aunque siempre existe 
una cota inferior trivial para cualquier algoritmo: al menos hay que leer los datos y 
luego escribirlos, de forma que ésa sería una primera cota inferior. Así, para 
ordenar n números una cota inferior sería n, y para multiplicar dos matrices de 
orden »n sería n”; sin embargo, los mejores algoritmos conocidos son de órdenes 
nlogn y n”” respectivamente. 


Propiedades de Q 


Veamos las propiedades de la cota inferior Q. La demostración de todas ellas se 
obtiene de forma simple aplicando la definición 1.2. 


1. Para cualquier función f'se tiene que fe Q(f). 

2. feQ(g) > Q(A) c Q(g). 

3. Q(/) =Q(g) S feQ(g) y g eQ(f). 

4. SifeQ(g) y g eQ(h) >f e Q(h). 

5. SifeQ(g) y fe Q(h) > f € Q(max(g,h)). 

6. Regla de la suma: Si f¡ e Q(g) y Ae Q(M > f, + fe Q(g + h). 
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7. Regla del producto: Si f¡ e Q(g) y f € Q(M) => f1:f2 € Q(g-h). 


8. Si existe iii ROD 
80 


a) Sik%+0 y k< oo entonces Q(f) = Q(g). 


b) Si k = 0 entonces g eQ(f), es decir, Q(g) < Q(f), pero sin embargo se 
verifica que f£ Q(g). 


= k, dependiendo de los valores que tome k obtenemos: 


De las propiedades anteriores se deduce que la relación —a, definida por f—a g 
si y sólo si Q(f) = Q(g), es una relación de equivalencia. Al igual que hacíamos 
para el caso de la cota superior O, siempre escogeremos el representante más 
sencillo para cada clase. Así los órdenes de complejidad (2 constante serán 
expresados por (1), los lineales por Q(n), etc. 


1.3.3 Orden Exacto. Notación O 


Como última cota asintótica, definiremos los conjuntos de funciones que crecen 
asintóticamente de la misma forma. 


Definición 1.3 
Sea f: N>][0,00). Se define el conjunto de funciones de orden O (Theta) de fcomo: 


0N=0070 
o, lo que es igual: 
OA = (g:N>10,00) | 3c,deR , c,d>0, 3nye N + cf) < e(m) < dfn) Vn > mo). 


Diremos que una función £ :N>[0,00) es de orden O de fsi £ e O(f). 


Intuitivamente, £ e O(f) indica que f está acotada tanto superior como 
inferiormente por múltiplos de f, es decir, que £ y fcrecen de la misma forma. 


Propiedades de O 


Veamos las propiedades de la cota exacta. La demostración de todas ellas se 
obtiene también de forma simple aplicando la definición 1.3 y las propiedades de O 
y Q. 


1. Para cualquier función f'se tiene que fe O(f). 

2. f£0(g8) > 0(/) = O(2). 

3. O) =0(28) SfE0(8) y g e€0(/. 

4. SifeO(g) y g€ € 0h) >f e O(h). 

5. Regla de la suma: Si f¡ € O(g) y £<0(M) > f¡ + f. € O(max(g,h)). 


10 TÉCNICAS DE DISEÑO DE ALGORITMOS 


6. Regla del producto: Si f¡ e€O(g) y A E 0(M) > ff € O(2»). 


7. Si existe lim Ja) 
>= g(n) 


a) Sik+0 y k< oo entonces O(f) = O(g). 
b) Sik=0 los órdenes exactos de f y g son distintos. 


= k, dependiendo de los valores que tome k obtenemos: 


1.3.4 Observaciones sobre las cotas asintóticas 


1. La utilización de las cotas asintóticas para comparar funciones de tiempo de 
ejecución se basa en la hipótesis de que son suficientes para decidir el mejor 
algoritmo, prescindiendo de las constantes de proporcionalidad. Sin embargo, 
esta hipótesis puede no ser cierta cuando el tamaño de la entrada es pequeño. 


2. Para un algoritmo dado se pueden obtener tres funciones que miden su tiempo 
de ejecución, que corresponden a sus casos mejor, medio y peor, y que 
denominaremos respectivamente 7, (1), T¡/2(n) y T,(n). Para cada una de ellas 
podemos dar tres cotas asintóticas de crecimiento, por lo que se obtiene un total 
de nueve cotas para el algoritmo. 


3. Para simplificar, dado un algoritmo diremos que su orden de complejidad es 
O(f) si su tiempo de ejecución para el peor caso es de orden O de f, es decir, 
T,(n) es de orden O(f). De forma análoga diremos que su orden de complejidad 
para el mejor caso es Q(g) si su tiempo de ejecución para el mejor caso es de 
orden (2 de g, es decir, 7, (1) es de orden Q(g). 


4. Por último, diremos que un algoritmo es de orden exacto O(f) si su tiempo de 
ejecución en el caso medio 7;/,(n) es de este orden. 


1.4 RESOLUCIÓN DE ECUACIONES EN RECURRENCIA 


En las secciones anteriores hemos descrito cómo determinar el tiempo de ejecución 
de un algoritmo a partir del cómputo de sus operaciones elementales (OE). En 
general, este cómputo se reduce a un mero ejercicio de cálculo. Sin embargo, para 
los algoritmos recursivos nos vamos a encontrar con una dificultad añadida, pues la 
función que establece su tiempo de ejecución viene dada por una ecuación en 
recurrencia, es decir, 7(n) = E(n), en donde en la expresión E aparece la propia 
función 7. 

Resolver tal tipo de ecuaciones consiste en encontrar una expresión no recursiva 
de 7, y por lo general no es una labor fácil. Lo que veremos en esta sección es 
cómo se pueden resolver algunos tipos concretos de ecuaciones en recurrencia, que 
son las que se dan con más frecuencia al estudiar el tiempo de ejecución de los 
algoritmos desarrollados según las técnicas aquí presentadas. 


1.4.1 Recurrencias homogéneas 


Son de la forma: 
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a T(n)+ a,T(n—1)+ a, T(n-2)+..+a,T(n-k)=0 


donde los coeficientes a, son números reales, y k es un número natural entre 1 y n. 
Para resolverlas vamos a buscar soluciones que sean combinaciones de funciones 
exponenciales de la forma: 


k 
T(n)=c,p,(n)r] +c,p,(mr3 +...+c,p, (mr =D cp (mr... 
i=1 


donde los valores c;, C»,...,Cn Y F1, F2, ...,, SON Números reales, y p¡(n),...,pi(n) son 
polinomios en n con coeficientes reales. Si bien es cierto que estas ecuaciones 
podrían tener soluciones más complejas que éstas, se conjetura que serían del 
mismo orden y por tanto no nos ocuparemos de ellas. 


Para resolverlas haremos el cambio x”" = T(n), con lo cual obtenemos la 
ecuación característica asociada: 


k k-1 k-2 
AX RH +a ox e a =0. 


Llamemos r;, r»,...,r a sus raices, ya sean reales o complejas. Dependiendo del 
orden de multiplicidad de tales raíces, pueden darse los dos siguientes casos. 


Caso 1: Raíces distintas 


Si todas las raíces de la ecuación característica son distintas, esto €s, 1,41; sii Xj, 
entonces la solución de la ecuación en recurrencia viene dada por la expresión: 


k 
n n n n 
T(n)=GK +0, +..+C 1H, = N cr, 
El 
donde los coeficientes c, se determinan a partir de las condiciones iniciales. 


Como ejemplo, veamos lo que ocurre para la ecuación en recurrencia definida 
para la sucesión de Fibonacci: 


T(n) = T(n-1) + T(n2), n>2 


con las condiciones iniciales 7(0) = 0, 7(1) = 1. Haciendo el cambio x= T(n) 
obtenemos su ecuación característica x” =x + 1, o lo que es igual, x*—x-—1=0, 
cuyas raíces son: 


14 4/5 1-45 
MA 
Z Z 
y por tanto 


e 


5 
T(n)=C, > +0, 3 
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Para calcular las constantes c, y cz necesitamos utilizar las condiciones iniciales 
de la ecuación original, obteniendo: 


mao Miedo 


T(0)=C, 3 HS 0. =c,+c,=0 , 
1 DC =C, ==. 
1+4/5 1/57. W/5 
TGD =, > +C, EN =1 


Sustituyendo entonces en la ecuación anterior, obtenemos 


E a EI : 
45 > EN > e O(p”). 


T(n)= 


Caso 2: Raices con multiplicidad mayor que 1 


Supongamos que alguna de las raíces (p.e. r¡) tiene multiplicidad m>1. Entonces la 
ecuación característica puede ser escrita en la forma 


m 
(—r Y" 742)..(0 — Fien+1) 
en cuyo caso la solución de la ecuación en recurrencia viene dada por la expresión: 


m k 


a il,,n n 
T()=) nr + Y ms 
¡=l 


i=m+1 


donde los coeficientes c,se determinan a partir de las condiciones iniciales. 
Veamos un ejemplo en el que la ecuación en recurrencia es: 


T(n) = 5T(n—-1) - 8T(n2) + 4T(n3), n22 
con las condiciones iniciales 7(k) = k para k=0, 1,2. La ecuación característica 
que se obtiene es x? — 5x? + 8x—-4=0, o lo que es igual (x-2)(x-1) =0 y por tanto, 
T(n) = c12" + can2” + c3l”. 
De las condiciones iniciales obtenemos c, = 2, c, =-1/2 y cz = 2, por lo que 
T(n) =2" — 12" -2 e0(n2”). 


Este caso puede ser generalizado de la siguiente forma. Si rj,r>,...,r son las 
raíces de la ecuación característica de una ecuación en recurrencia homogénea, 
cada una de multiplicidad m,, esto es, si la ecuación característica puede expresarse 
como: 


(ar) (ar). (xr) =0, 


entonces la solución a la ecuación en recurrencia viene dada por la expresión: 
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mm ma Mp 
¡1 ¡1 ¡1 
T(n) = ) CATE +) CN EH + ) EM 
¡=1 i=l ¡=1 


1.4.2 Recurrencias no homogéneas 


Consideremos una ecuación de la forma: 
a,T(m+a T(n—-1)+..+a,T(n—-k)=b"p(n) 


donde los coeficientes a; y b son números reales, y p(n) es un polinomio en n de 
grado d. Una primera idea para resolver la ecuación es manipularla para convertirla 
en homogénea, como muestra el siguiente ejemplo. 


Sea la ecuación 7(n) — 27(n-1) = 3” para n > 2, con las condiciones iniciales 
T(0) =0 y 7(1) = 1. En este caso b= 3 y p(n) = 1, polinomio en n de grado 0. 


Podemos escribir la ecuación de dos formas distintas. En primer lugar, para n+1 
tenemos que 


T(n+1)-2T(n)= 3"". 
Pero si multiplicamos por 3 la ecuación original obtenemos: 
3T(n) -6T(n-1)=3"" 
Restando ambas ecuaciones, conseguimos 
T(n+1) -5T(n) + 67(n-1) =0, 


que resulta ser una ecuación homogénea cuya solución, aplicando lo visto 
anteriormente, es 


T(n)= 3" - 2" eO(3”). 


Estos cambios son, en general, difíciles de ver. Afortunadamente, para este tipo 
de ecuaciones también existe una fórmula general para resolverlas, buscando sus 
soluciones entre las funciones que son combinaciones lineales de exponenciales, en 
donde se demuestra que la ecuación característica es de la forma: 


(ax +ax" +ax +..+9,0—-by" =0, 


lo que permite resolver el problema de forma similar a los casos anteriores. 


Como ejemplo, veamos cómo se resuelve la ecuación en recurrencia que plantea 
el algoritmo de las torres de Hanoi: 
T(n) =2T(n-1) +1. 
Su ecuación característica es entonces (x-2)-1)? = 0, y por tanto 
T(n) = c,2" + c21" +c3n1" e€6Q”). 


Generalizando este proceso, supongamos ahora una ecuación de la forma: 
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a,T(m+a T(n-1)+..+a,T(n—-k)=b/p,(nM)+b) p,(M)+...+b/; p (n) 


donde como en el caso anterior, los coeficientes a; y b, son números reales y p(n) 
son polinomios en n de grado d),. En este caso también existe una forma general de 
la solución, en donde se demuestra que la ecuación característica es: 


(araxira O a ba oy oa y? ab y" =0. 


Como ejemplo, supongamos la ecuación 
T(n) =2T(n-1) +n+2", n21, 


con la condición inicial 7(0) = 1. En este caso tenemos que bi = L pin) = 
b2=2 y pin) =1, por lo que su ecuación característica es (2-1) =0, lo ans 
da lugar a la expresión final de 7(m): 


T(n)=2-n +2 +2" eQ(n2”). 


1.4.3 Cambio de Variable 
Esta técnica se aplica cuando n es potencia de un número real a, esto es, n = a”. Sea 
por ejemplo, para el caso a = 2, la ecuación T(n) = 47(n/2) + n, donde n es una 
potencia de 2 (n > 3), 7(1) = 1, y 72) = 6. 
Si n= 2% podemos escribir la ecuación como: 

TP =4T104%) +2* 
Haciendo el cambio de variable £; = 7(2*) obtenemos la ecuación 

t¿= ty + 2% 


que corresponde a una de las ecuaciones estudiadas anteriormente, cuya solución 
viene dada por la expresión 


1 = 012%) + 0,2% 

Deshaciendo el cambio que realizamos al principio obtenemos que 
T(n) = ci” + cm. 

Calculando entonces las constantes a partir de las condiciones iniciales: 


T(n)=2n?-—n e 0(n). 


1.4.4 Recurrencias No Lineales 


En este caso, la ecuación que relaciona 7(n) con el resto de los términos no es 
lineal. Para resolverla intentaremos convertirla en una ecuación lineal como las que 
hemos estudiado hasta el momento. 


Por ejemplo, sea la ecuación T(n) = nT “(n/2) para n potencia de 2,n > 1, con la 
condición inicial 7(1) = 1/3. Llamando 1¿= 7(2*), la ecuación queda como 
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t, = T(2%) e A = ee 


que no corresponde a ninguno de los tipos estudiados. Necesitamos hacer un 
cambio más para transformar la ecuación. Tomando logaritmos a ambos lados y 
haciendo el cambio u,= log t, obtenemos 


U;— 241 = K, 


ecuación en recurrencia no homogénea cuya ecuación característica asociada es 
(6-2)0—1)= 0. Por tanto, 


UL c¡2* + Ca + C3k. 
Necesitamos ahora deshacer los cambios hechos. Primero u;= log f; 
_D 2% +c, +cz3k 
t,=2 
y después 1¿= T(2*). En consecuencia 


logn 
T(n) =3 yr +0) +C% : 


Para calcular las constantes necesitamos las condiciones iniciales. Como sólo 
disponemos de una y tenemos tres incógnitas, usamos la ecuación en recurrencia 
original para obtener las restantes: 


T(2)=2TU1)=2/9. 
T(4) =4T*(0) = 16/81. 


Con esto llegamos a que c¡ = log(4/3) = 2 — log3, cz = 2, c3 = —1 y por 
consiguiente: 


92 
4n3" 


Tn) = 


16 
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1.5 PROBLEMAS PROPUESTOS 


1.1. 


1.2. 


1.3. 


1.4. 


De las siguientes afirmaciones, indicar cuales son ciertas y cuales no: 


(1)  nre0(n) (ix) meQ() 

(id) ne0(n) () reQ(í) 

(iii) 22000”) (a) 200”) 

(iv) (n+1)!e O(n!) (ii) (n+D)leQ(n!) 

(y  Afme0m => 2"e0(2”) (xiid) An) Q(m) > 2 Me Q(2”) 
(vi) 3"e0(2”) (xiv) 3"eQ(Q”) 

(vii) logn eO(n!”) (xv) logn e Q(n'?) 

(viii) n'?O(logn) (xvi) n'eQílogn) 


Sea a una constante real, 0<a<1. Usar las relaciones € y = para ordenar los 
. +a 


órdenes de complejidad de las siguientes funciones: nlogn, nlogn, n*, n *, 
(1+a), (1*+8n+logóny", n"/logn, 2”. 


La siguiente ecuación recurrente representa un caso típico de un algoritmo 


recursivo: 
k . 
cn sil<n<b 
T(n)= O 
aT(n-b)j+cn" sin>b 


donde a,c,k son números reales, n,b son números naturales, y a>0, c>0, 
k20. En general, la constante a representa el número de llamadas recursivas 
que se realizan para un problema de tamaño n en cada ejecución del 
algoritmo; n-b es el tamaño de los subproblemas generados; y cn 
representa el coste de las instrucciones del algoritmo que no son llamadas 
recursivas. 


On)  sia<l 


ci 


Demostrar que T(n)e ¿ O(n sia=1 


O(a"“?) sia>1 


La siguiente ecuación recurrente representa un caso típico de Divide y 
Vencerás: 


k 


cn sil<n<b 
T(n)= Pa 
aT(n/b)+cn" sin2b 


donde a,c,k son números reales, n,b son números naturales, y a>0, c>0, 
k20, b>1. La expresión cn” representa en general el coste de descomponer 
el problema inicial en a subproblemas y el de componer las soluciones para 
producir la solución del problema original. 
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On) sia<b* 
Demostrar que T(n)e¿0O(n' logn) sia=b* 
One) sia>b* 


1.5. Supongamos que disponemos de la siguiente definición de tipo: 


CONST n= ...; 
TYPE vector = ARRAY [1..n] OF INTEGER; 


Consideramos entonces los procedimientos y funciones siguientes: 
PROCEDURE Algoritmo1(VAR a:vector)'; 


VAR i,j:CARDINAL; 
temp: INTEGER; 


BEGIN 
FOR i:=1 TO n-1 DO (* 1 x) 
FOR j:=n TO ¿+1 BY -1 DO (* 2 x) 
IF a[j-11>a[j] THEN (* 3 x) 
temp:=a[j-11; (* 4 x) 
a[j-11:=a[j]; (* 5 x) 
a[jl:=temp (* 6 x) 
END (* 7 x) 
END (* 8 x) 
END (* 9 *) 


END Algoritmol1; 


PROCEDURE Algoritmo2(VAR a:vector;c: INTEGER) : CARDINAL; 
VAR inf,sup,i:CARDINAL; 


BEGIN 
inf:=1; sup:=n; (x* 1 *) 
WHILE (sup>=inf) DO (* 2 x) 
i:=(inf+sup) DIV 2; (* 3 x) 
IF a[lil=c THEN RETURN i (* 4 x) 
ELSIF c<a[il] THEN sup:=i-1 (* 5 x) 
ELSE inf:=i+1 (* 6 x) 
END (* 7 x) 
END; (x* 8 x) 
RETURN 0; (x* 9 x) 


END Algoritmo2; 
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PROCEDURE Euclides (m,n: CARDINAL) : CARDINAL; 
VAR temp: CARDINAL; 


BEGIN 
WHILE m>0 DO (* 1 x) 
temp:=m; (x* 2 *) 
m:=n MOD m; (* 3 x) 
n:=temp (* 4 1) 
END; (* 5 x) 
RETURN n (x* 6 x) 
END Euclides; 
PROCEDURE Misterio(n:CARDINAL) ; 
VAR 1,j,k,s: INTEGER; 
BEGIN 
s:=0; (* 1 x) 
FOR i:=1 TO n-1 DO (* 2 x*) 
FOR j:=i+1 TO n DO (x* 3 *) 
FOR k:=1 TO j DO (* 4 x) 
s:=s+2 (* 5 x) 
END (x* 6 x) 
END (xk 7 x) 
END (* 8 x) 


END Misterio; 


a) Calcular sus tiempos de ejecución en el mejor, peor, y caso medio. 
b) Dar cotas asintóticas O, (Q y O para las funciones anteriores. 


1.6. Demostrar las siguientes inclusiones estrictas: O(1) € O(logn) € O(n) < 
O(nlogn) < O(n”) < O(1*) < O(nÉ) < O(2”) c O(n!). 


1.7. a) Demostrar que fe O(g) S g e Q(f). 
b) Dar un ejemplo de funciones f y g tales que fe O(g) pero que f£ Q(g). 
c) Demostrar que Va,b>1 se tiene que log,ne O(log,n). 


1.8.  Considérense las siguientes funciones de »: 
Am =n, f(n)=n* + 1000n; 
0 n,  sinimpar Li n, sins100 
n)= s n)= 
7 n?, sinpar i n”, sin>100 


Para cada posible valor de ¡,j indicar si fe O(f)) y si fe Q(f). 
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1.9. 


1.10. 


Resolver las siguientes ecuaciones y dar su orden de complejidad: 


a) T(M=3T(n-1)+4T(n2) sin>1; 7(0)50; T(D=1. 

b) T(n)=2T(n-1)An+5)y3" si 12>0; 7(0)-0. 

c) T(M-=4T(1/2 yen? si 1>4, n potencia de 2; 7(1)=1; 72)-8. 
d) T(M=2T(1/2)+nlogn si n>1l, n potencia de 2. 

e) T(1M)+3T(n/2)+5n+3 si n>1l, n potencia de 2. 

D Tmy2TM/2)Hogn si n>1l, n potencia de 2. 

2) Tmy-2T(1 2) +Hogn con n=22, TOF1. 

h) T(m)=5T(1/2)+(mlogny si n>1,n potencia de 2; 7(1)=1. 


1) T)=T(n-1)92T(n2)2T(n-3) si n>2; T(n)=9n"-15n+106 si n=0,1,2. 
Y Tm-6/2370/2-1/2370/41/n) si n>2; T(D)=1; 70)+3/2. 

kx) TM-2T(n/4)+n'P si 1>4, n potencia de 4. 

1) T(m)=4T(0/3)+n? s1 1>3, n potencia de 3. 


Suponiendo que 7,8 O(f) y que T,e O(f), indicar cuáles de las siguientes 
afirmaciones son ciertas: 


a) T,+T,€0(/. 
b) 7, -T,E€0(f. 
c) 7,/T,€0(1). 
d) T,€0(T7). 


Encontrar dos funciones (1) y g(n) tales que f£ O(g) y g £0(f). 
Demostrar que para cualquier constante k se verifica que log'n e O(n). 


Consideremos los siguientes procedimientos y funciones sobre árboles. 
Calcular sus tiempos de ejecución y sus órdenes de complejidad. 


PROCEDURE Inorden(t:arbol); (* recorrido en inorden de t x*) 
BEGIN 
IF NOT Esvacio(t) THEN Ge 1 
*) 
Inorden(Izq(t)); (x* 2 x) 
Opera(Raiz(t)); (x* 3 *) 
Inorden (Der (t)) (x* 4 x) 
END; (x* 5 x) 


END Inorden; 
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PROCEDURE Altura(t:arbol):CARDINAL; (* altura de t *) 


BEGIN 
IF Esvacio(t) THEN (* 1 x) 
RETURN O (x* 2 x) 
ELSE (x* 3 x) 
RETURN 1+Max2(Altura(Izq(t)),Altura(Der(t))) (x 4 x) 
END (x* 5 x) 
END Altura; 


PROCEDURE Mezcla(t1,t2:arbol):arbol; 

(* devuelve un arbol binario de busqueda con los elementos de 
los dos arboles binarios de busqueda ti y t2. La funcion Ins 
inserta un elemento en un arbol binario de busqueda *) 


BEGIN 
IF Esvacio(t1) THEN (x* 1 x) 
RETURN t2 (* 2 *) 
ELSIF Esvacio(t2) THEN (* 3 x) 
RETURN t1 (* 4 x) 
ELSE (* 5 *) 
RETURN Mezcla(Mezcla(Ins(t1,Raiz(t2)),Izq(t2)), 
Der (t2)) (x* 6 x) 
END (* 7 x) 
END Mezcla; 


Supondremos que las operaciones básicas del tipo abstracto de datos arbol 
(Raiz, Izq, Der, Esvacio) son O(1), así como las operaciones Opera 
(que no es relevante lo que hace) y Max2 (que calcula el máximo de dos 
números). Por otro lado, supondremos que la complejidad de la función Ins 
es O(logn). 


1.14. Ordenar las siguientes funciones de acuerdo a su velocidad de crecimiento: 
n, Jn, logn, loglogn, log?n, n/logn, Án log?n, (1/3), (3/2)”, 17, n?. 


1(€ A 

1.15. Resolver la ecuación 7 (n) = ($ 10) + cn, siendo 7 (0) =0. 
AX i=o 

1.16. Consideremos las siguientes funciones: 


CONST n= ...; 
TYPE vector = ARRAY[1..n] OF INTEGER; 
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PROCEDURE BuscBin(VAR a:vector; 
prim,ult:CARDINAL;x: INTEGER) : BOOLEAN; 
VAR mitad:CARDINAL; 


BEGIN 
IF (prim>=ult) THEN RETURN a[ult]=x (* 1 x) 
ELSE (* 2 x) 
mitad:=(prim+ult)DIV 2; (* 3 x) 
IF x=a[mitad] THEN RETURN TRUE (* 4 x) 
ELSIF (x<a[mitadl) THEN (* 5 x) 
RETURN BuscBin(a,prim,mitad-1,x) (* 6 x) 
ELSE (* 7 x) 
RETURN BuscBin(a,mitad+1,ult,x) (x 8 x) 
END (* 9 x) 
END (* 10 x) 
END BuscBin; 
PROCEDURE Sumadigitos (num: CARDINAL) : CARDINAL; 
BEGIN 
IF num<1í0 THEN RETURN num (x* 1 x) 
ELSE RETURN (num MOD 10)+Sumadigitos(num DIV 10) (x* 2 x*) 
END (x* 3 x) 


END Sumadigitos; 


a) Calcular sus tiempos de ejecución y sus órdenes de complejidad. 

b) Modificar los algoritmos eliminando la recursión. 

c) Calcular la complejidad de los algoritmos modificados y justificar para 
qué casos es más conveniente usar uno u otro. 


1.17. Consideremos la siguiente función: 


PROCEDURE Raro(VAR a:vector;prim,ult:CARDINAL) : INTEGER; 
VAR mitad,terc:CARDINAL; 


BEGIN 
IF (prim>=ult) THEN RETURN a[ult] END; 
mitad:=(prim+ult)DIV 2; (* posicion central *) 
terc :=(ult-prim)DIV 3; (* num. elementos DIV 3 *) 
RETURN a[mitad]+Raro(a,prim,prim+terc)+Raro(a,ult-terc,ult) 
END Raro; 


a) Calcular el tiempo de ejecución de la llamada a la función Raro(a,1,n), 
suponiendo que n es potencia de 3 . 


b) Dar una cota de complejidad para dicho tiempo de ejecución. 
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1.6 SOLUCIÓN A LOS PROBLEMAS PROPUESTOS 


Antes de comenzar con la resolución de los problemas es necesario hacer una 
aclaración sobre la notación utilizada para las funciones logarítmicas. A partir de 
ahora y a menos que se exprese explícitamente otra base, la función “log” hará 
referencia a logaritmos en base dos. 


Solución al Problema 1.1 (9/8) 
(1)  n?e0(n) es cierto pues lim, (nn) =0. 
(ii) nie O(n) es falso pues lim, _.. (1*/n7) =0. 


(iii) 2""0(2”) es cierto pues lim,_,..(2"""/2")=2. 


miren 


(iv) (n+1)!eO(n!) es falso pues lim, _.. (n!/(n+1)!) =0. 


(y) fime0O(n) > 2"e0(2”) es falso. Por ejemplo, sea (nm) = 3n; claramente 
Am) O(n) pero sin embargo lim, _._ (2”/2%) = 0, con lo cual 2” O(2”). De 


forma más general, resulta ser falso para cualquier función lineal de la forma 
fin) = an con > 1, y cierto para fín) = Pn con P< 1. 


n>0 


(vi) 3"£0(2”) es falso pues lim,_...(2"/3”) =0. 


(vii) logn eO(n!?) es cierto pues lim, (logn/n'?) =0. 


n>0o 


(viii) n'?e O(logn) es falso pues lim, _._. (logn/n'”?) =0. 


(ix) nr eQ(n) es falso pues lim,_... (n/m) =0. 
(x) ne Q(1) es cierto pues lim, (nn?) =0, 


(xi) 2"*eQ(”) es cierto pues lim __(2""/2”) =2. 


n—o0 


(xii) (n+1)!le Q(n!) es cierto pues lim, __ (n!/(n+1)!) =0. 


n>o 

(xiii) Ane Q(m) => 2" Q(”) es falso. Por ejemplo, sea (nm) = (1/2)m; claramente 
An)e O(n) pero sin embargo lim, _..(20?”/2”) =0, con lo cual 22 Q(2”). 
De forma más general, resulta ser falso para cualquier función f(n) = On con 
0a< 1, y cierto para in) = Bn con P> 1. 


(xiv) 3"eQ(”) es cierto pues lim, ,., (2"/3")=0. 
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1/2 


(xv) logn e Q(n!”) es falso pues lim, (logn/n'?) =0. 


n—>o0 


1/2 


(xvi) n'“e Q(logn) es cierto pues lim, _.. (logn/n'?) =0. 


Solución al Problema 1.2 (O) 
e Respecto al orden de complejidad O tenemos que: 


O(nlogn) < O(n'**) c O(n 1logn) < O(n logn) < O(n*) = O((n+8n+log ny) 
cO((1+a)) c 02”). 
Puesto que todas las funciones son continuas, para comprobar que O(f)CO(g), 
basta ver que lim,_.. (An)g(n)) = 0, y para comprobar que O(f) = O(g), basta 
ver que lim, _... (Am)/g(n)) es finito y distinto de 0. 


e Por otro lado, respecto al orden de complejidad Q, obtenemos que: 
Qí(nlogn) > Q(n'***”) > Q(1logn) > Q(nélogn) > Q(a?) = Q((n*+8n+log'ny*) > 
Q(U1+a)”) > QQ”) 
Para comprobar que Q(f) C Q(g), basta ver que lim, _...(g(n)/fKn)) = 0, y para 
comprobar que Q(f) = Q(g), basta ver que lim, (An)/g(n)) es finito y distinto 


de 0 puesto que al ser las funciones continuas tenemos garantizada la existencia 
de los límites. 


e Y en lo relativo al orden de complejidad O, al definirse como la intersección de 
los órdenes O y Q, sólo tenemos asegurado que: 


O(n*) = O((1+8n+log ny), 


siendo los órdenes O del resto de las funciones conjuntos no comparables. 


Solución al Problema 1.3 (ar) 


7 “7 . k a 
La ecuación dada puede ser también escrita como T(n) — aT(n-b) = cn”, ecuación 
en recurrencia no homogénea cuya ecuación característica es: 


La) iy =0. 


e Para estudiar las raíces de esa ecuación, vamos a suponer primero que a + 1. En 
este caso, la ecuación tiene una raíz de multiplicidad k+1 (el 1), y b raíces 
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y e e a F ., 
distintas ry,7,...,rp (las b raíces b-ésimas de a' ). Entonces la solución de la 
ecuación en recurrencia es de la forma: 


T(M=c,V" +e,n1" +01 +..+0,¿N 1 +d 8 +d,03 +..+d,rj = 


k+1 


12) 
¡1 
Hen Yarz 
i=1 i=1 


siendo c; y d; coeficientes reales. 


— Sia<l, las raíces b-ésimas de a (esto es, las r;) son menores en módulo que 
1, con lo cual el segundo sumatorio tiende a cero cuando n tiende a co, y en 


y . Tn , ie 
consecuencia 7(1) e O(n*) pues lim sito) = Cy, es finito y distinto de cero. 
n>o y 
Para ver que efectivamente Cy, es distinto de cero independientemente de las 
condiciones iniciales, sustituimos esta expresión de 7(n) en la ecuación 
original, 7(n) — aT(n-b) = cn”, y por tanto: 


k+1 b k+l b 


NO -a Nano + de om. 
i=l i=l i=l i=l 


a m k 
Igualando ahora los coeficientes que acompañan a n* obtenemos que 
Ck+1 — ACh+1 = C, O lo que es igual, (1-a)c;+¡ = c. Ahora bien, como sabemos 
que a< 1 y c> 0, entonces cCz+/ no puede ser cero. 


— Sia> l, las funciones del segundo sumatorio son exponenciales, mientras 
que las primeras se mantienen dentro de un orden polinomial, por lo que en 
este caso el orden de complejidad del algoritmo es exponencial. Ahora bien, 
como todas las raíces b-ésimas de a tienen el mismo módulo, todas crecen de 
la misma forma y por tanto todas son del mismo orden de complejidad, 
obteniendo que 


00) =0(7)=...=0(7). 


Como lim EU) 


n—0 A 


= dl, es distinto de cero y finito, podemos concluir que: 


T(n)e O(1")= o (Va) ) =O(ar o), 


Hemos supuesto que d, + 0. Esto no tiene por qué ser necesariamente cierto 
para todas las condiciones iniciales, aunque sin embargo sí es cierto que al 
menos uno de los coeficientes d, ha de ser distinto de cero. Basta tomar ese 
sumando para demostrar lo anterior. 


* Recordemos que dados dos números reales a y b, la solución de la ecuación x”— a =0 
j É ae 1/b 2mik) 
tiene hb raíces distintas, que pueden ser expresadas como a “e”””, para (=0,1,2,...,n—1. 
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e Supongamos ahora que a = 1. En este caso la multiplicidad de la raíz 1 es k+2, 
con lo cual 


k+1 


Tm =c,P +on"+cnó 1 +. +eÓ ni V 4d +... +d,r] = 


(Len )e( Lar 


Pero las raíces r,,r3,...,rp son todas de módulo 1 (obsérvese que r¡=1), y por 
tanto el segundo sumando de 7(n) es de complejidad O(1). 


Así, el crecimiento de 7(n) coincide con el ¡del primer sumando, que es un 
polinomio de grado k+1 con lo cual T(m)e O(n**?). 


Solución al Problema 1.4 (6) 
Haciendo el cambio n = b”, o lo que es igual, m = log¿n, obtenemos que 
T(b”) =aT(b"") + cb”. 
Llamando t,,= 7(b”), la ecuación queda como 
Em Ata > abty”, 


ecuación en recurrencia no homogénea con ecuación característica (x-aXx—b") = 0. 


k 
Para resolver esta ecuación supongamos primero que a = b”. Entonces, la 
ecuación característica es (x—b')?= 0 y por tanto 


tn =C¡b"” + eomb*”. 
Necesitamos ahora deshacer los cambios hechos. Primero t,,= 7(b”) con lo que 
T(b”) = cb" + comb'” = (c+ camyb"”, 
y después n = b”, obteniendo finalmente que 
T(n) = (c, + calog¿n)n! e O(n'logn).* 


Supongamos ahora el caso contrario, a + b*. Entonces la ecuación característica 
tiene dos raíces distintas, y por tanto 


P=arraco", 
Necesitamos deshacer los cambios hechos. Primero f,, = 7(b”), con lo que 
T(b”) = cia” + cab?”, 
y después n = b”, obteniendo finalmente que 


peeva 


1 k 
Tm=c ¡a "cm =cau" +0, 


* Obsérvese que se hace uso de que log¿ne O(logn), lo que se demuestra en el problema 
1.7. 
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y . Ñ 7 . 1 . 
En consecuencia, si log,a > k (si y sólo si a> b*) entonces T(n) e An” *). Si 
no, es decir, a < b”, entonces T(n) e O(m5). 


Solución al Problema 1.5 


Procedimiento Algoritmol (O) 


a) Para obtener el tiempo de ejecución, calcularemos primero el número de 
operaciones elementales (OE) que se realizan: 


— En la línea (1) se ejecutan 3 OE (una asignación, una resta y una comparación) 
en cada una de las iteraciones del bucle más otras 3 al final, cuando se efectúa la 
salida del FOR. 


— Igual ocurre con la línea (2), también con 3 OE (una asignación, una suma y una 
comparación) por iteración, más otras 3 al final del bucle. 


— En la línea (3) se efectúa una condición, con un total de 4 OE (una diferencia, 
dos accesos a un vector, y una comparación). 


— Las líneas (4) a (6) sólo se ejecutan si se cumple la condición de la línea (3), y 
realizan un total de 9 OE: 3, 4 y 2 respectivamente. 


Con esto: 


En el caso mejor para el algoritmo la condición de la línea (3) será siempre 
falsa, y no se ejecutarán nunca las líneas (4), (5) y (6). Así, el bucle más interno 
realizará (n—1) iteraciones, cada una de ellas con 4 OE (línea 3), más las 3 OE de 
la línea (2). Por tanto, el bucle más interno realiza un total de 


Y (443) |+3=7] Y 1|+3=7(n-1)+3 
j=i+1 j=i+1 
OE, siendo el 3 adicional por la condición de salida del bucle. 


A su vez, el bucle externo repetirá esas 7(n-1)+3 OE en cada iteración, lo 
que hace que el número de OE que se realizan en el algoritmo sea: 


i=1l 


T(n) = (E 0=0+3)03)+3 20 203 


e En el caso peor, la condición de la línea (3) será siempre verdadera, y las líneas 
(4), (5) y (6) se ejecutarán en todas las iteraciones. Por tanto, el bucle más 
interno realiza 


[Euros Jes=tón=0. 


j=i+1 


OE. El bucle externo realiza aquí el mismo número de iteraciones que en el caso 
anterior, por lo que el número de OE en este caso es: 
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T(n) = (Eom—0 +3) ese -2n-3. 


i=l 


e Enel caso medio, la condición de la línea (3) será verdadera con probabilidad 
1/2. Así, las líneas (4), (5) y (6) se ejecutarán en la mitad de las iteraciones del 
bucle más interno, y por tanto realiza 


[Ear jr dao 


j=i+l 


OE. El bucle externo realiza aquí el mismo número de iteraciones que en el caso 
anterior, por lo que el número de OE en este caso es: 


n—-1 
T(n) = (E Po=0+3)+3)e3 ho +2n-3, 


i=1 


b) Como los tiempos de ejecución en los tres casos son polinomios de grado 2, la 
complejidad del algoritmo es cuadrática, independientemente de qué caso se trate. 


Obsérvese cómo hemos analizado el tiempo de ejecución del algoritmo sólo en 
función de su código y no respecto a lo que hace, puesto que en muchos casos esto 
nos llevaría a conclusiones erróneas. Debe ser a posteriori cuando se analice el 
objetivo para el que fue diseñado el algoritmo. 


En el caso que nos ocupa, un examen más detallado del código del 
procedimiento nos muestra que el algoritmo está diseñado para ordenar de forma 
creciente el vector que se le pasa como parámetro, siguiendo el método de la 
Burbuja. Lo que acabamos de ver es que sus casos mejor, peor y medio se 
producen respectivamente cuando el vector está inicialmente ordenado de forma 
creciente, decreciente y aleatoria. 


Función A/goritmo2 (és) 


a) Para calcular el tiempo de ejecución, calcularemos primero el número de 
operaciones elementales (OE) que se realizan: 


— En la línea (1) se ejecutan 2 OE (dos asignaciones). 

— En la línea (2) se efectúa la condición del bucle, que supone 1 OE (la 
comparación). 

— Las líneas (3) a (6) componen el cuerpo del bucle, y contabilizan 3, 2+1, 2+2 y 
2 OE respectivamente. Es importante hacer notar que el bucle también puede 
finalizar si se verifica la condición de la línea (4). 


— Por último, la línea (9) supone 1 OE. A ella se llega cuando la condición del 
bucle WHILE deja de ser cierta. 


Con esto: 
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c) 
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En el caso mejor se efectuarán solamente la líneas (1), (2), (3) y (4). En 
consecuencia, 7(n) = 2+1+3+3 =9, 


En el caso peor se efectúa la línea (1), y después se repite el bucle hasta que su 
condición sea falsa, acabando la función al ejecutarse la línea (9). Cada 
iteración del bucle está compuesta por las líneas (2) a (8), junto con una 
ejecución adicional de la línea (2) que es la que ocasiona la salida del bucle. En 
cada iteración se reducen a la mitad los elementos a considerar, por lo que el 
bucle se repite logn veces. Por tanto, 


log n 

1-24 [[Ha+r1922+2)+1)e1=10108n+4 
i=l 

En el caso medio, necesitamos calcular el número medio de veces que se repite 

el bucle, y para esto veamos cuántas veces puede repetirse, y qué probabilidad 

tiene cada una de suceder. 


Por un lado, el bucle puede repetirse desde una vez hasta logn veces, puesto 
que en cada iteración se divide por dos el número de elementos considerados. Si 
se repitiese una sola vez, es que el elemento ocuparía la posición n/2, lo que 
ocurre con una probabilidad 1/(1+1). Si el bucle se repitiese dos veces es que el 
elemento ocuparía alguna de las posiciones n/4 ó 3n/4, lo cual ocurre con 
probabilidad 1/(1+1)+1/(1+1)=2/(n+1). En general, si se repitiese ¡ veces es que 
el elemento ocuparía alguna de las posiciones nk/2', con k impar y 1<k<2". 

Es decir, el bucle se repite ¡ veces con probabilidad 2*/(n+1). Por tanto, el 
número medio de veces que se repite el ciclo vendrá dado por la expresión: 


E. 27  nlogn—n+1 
1 1 
la n+l n+1 


Con esto, la función ejecuta la línea (1) y después el bucle se repite ese 
número medio de veces, saliendo por la instrucción RETURN en la línea (4). 
Por consiguiente, 
nlogn—n+1 

n+1 


O DE 
n+ 


En el caso mejor el tiempo de ejecución es una constante. Para los casos peor y 
medio, la complejidad resultante es de orden O(logn) puesto que 


es una constante finita y distinta de cero en ambos casos (10 y 8 
respectivamente). 
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Función Euclides (8) 


a) En este caso el análisis del tiempo de ejecución y la complejidad de la función 
sigue un proceso distinto al estudiado en los casos anteriores. 


Lo primero es resaltar algunas características del algoritmo, siguiendo una línea 
de razonamiento similar a la de [BRA97]: 


[1] Para cualquier par de enteros no negativos m y n tales que n>m, se verifica que 
n MOD m< n/2. Veámoslo: 


a) Sim > n/2 entonces 1< n/m < 2 y por tanto n DIV m = 1, lo que implica 
que n MOD m=n— mín DIV m) =n-m<n-—n/2=n/2. 


b) Por otro lado, si m < n/2 entonces n MOD m <m< n/2. 


[2] Podemos suponer sin pérdida de generalidad que n > m. Si no, la primera 
iteración del bucle intercambia n con m ya que n MOD m = nm cuando n < m. 
Además, la condición n > m se conserva siempre (es decir, es un invariante del 
bucle) pues n MOD m nunca es mayor que m. 


[3] El cuerpo del bucle efectúa 4 OE, con lo cual el tiempo del algoritmo es del 
orden exacto del número de iteraciones que realiza el bucle. Por consiguiente, 
para determinar la complejidad del algoritmo es suficiente acotar este número. 


[4] Una propiedad curiosa de este algoritmo es que no se produce un avance 
notable con cada iteración del bucle, sino que esto ocurre cada dos iteraciones. 
Consideremos lo que les ocurre a m y n cuando el ciclo se repite dos veces, 
suponiendo que no acaba antes. Sean mo y no los valores originales de los 
parámetros, que podemos suponer ny > my por [2]. Después de la primera 
iteración, m vale ny MOD my. Después de la segunda iteración, n toma ese 
valor, y por tanto ya es menor que ny/2 (por [1]). En consecuencia, n vale 
menos de la mitad de lo que valía tras dos iteraciones del bucle. Como se 
sigue manteniendo que n > m, el mismo razonamiento se puede repetir para las 
siguientes dos iteraciones, y así sucesivamente. 


El hecho de que n valga menos de la mitad cada dos iteraciones del bucle es el 
que nos permite intuir que el bucle se va a repetir del orden de 2logn veces. Vamos 
a demostrar esto formalmente. 


Para ello, vamos a tratar el bucle como si fuera un algoritmo recursivo. Sea 7(/) 
el número máximo de veces que se repite el bucle para valores iniciales m y n 
cuando m < n < 1. En este caso / representa el tamaño de la entrada. 


— Sin<2el bucle no se repite (si m= 0) o se hace una sola vez (si mes 1 6 2). 
— Sin>2ym=1 0 bien m divide a n, el bucle se repite una sola vez. 
— En otro caso (n > 2 y m no divide a n) el bucle se ejecuta dos veces, y por lo 


visto en [4], n vale a lo sumo la mitad de lo que valía inicialmente. En 
consecuencia n < (1/2), y además m se sigue manteniendo por debajo de n. 


Esto nos lleva a la ecuación en recurrencia 7(/) <2 + 7(1/2) s11>2, (1) <1 si/< 
2, lo que implica que el algoritmo de Euclides es de complejidad logarítmica 
respecto al tamaño de la entrada (/). 
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Nos preguntaremos la razón de usar 7(/) para acotar el número de iteraciones 
que realiza el algoritmo en vez de definir T directamente como una función de n, el 
mayor de los dos operandos, lo cual sería mucho más intuitivo. 


El problema es que si definimos 7(nm) como el número de iteraciones que realiza 
el algoritmo para los valores m < n, no podríamos concluir que 7(n) € 2 + T(n/2) 
del hecho de que n valga la mitad de su valor tras cada dos iteraciones del bucle. 


Por ejemplo, para Euclides(8,13), obtenemos que 7(13) = 5 en el peor caso, 
mientras que 7(13/2) = 7(6) = 2. Esto ocurre porque tras dos iteraciones del bucle n 
no vale 6, sino 5 (y m= 3), y con esto sí es cierto que 7(13) < 2 + T(5) ya que 7(5) 
=3, 

La raíz de este problema es que esta nueva definición más intuitiva de 7 no 
lleva a una función monótona no decreciente (7(5) > 7(6)) y por tanto la existencia 
de algún n* < n/2 tal que T(n) < 2 + T(n”) no implica necesariamente que 
T(n) <2 + T(n/2). 

En vez de esto, solamente podríamos afirmar que 7(1)£2+max(T(n”)|n"<n2 y, 
que es una ecuación en recurrencia bastante difícil de resolver. Esa es la razón de 
que escogiésemos nuestra función 7 de forma que fuera no decreciente y que 
expresara una cota superior del número de iteraciones. 


Para acabar, es interesante hacer notar una característica curiosa de este 
algoritmo: se demuestra que su caso peor ocurre cuando m y n son dos términos 
consecutivos de la sucesión de Fibonacci. 


b) T(D)e O(log/) como se deduce de la ecuación en recurencia que define el tiempo 
de ejecución del algoritmo. 


Procedimiento Misterio (O) 


a) En este caso son tres bucles anidados los que se ejecutan, independientemente de 
los valores de la entrada, es decir, no existe peor, medio o mejor caso, sino un 
Único caso. 

Para calcular el tiempo de ejecución, veamos el número de operaciones 
elementales (OE) que se realizan: 


— En la línea (1) se ejecuta 1 OE (una asignación). 


— En la línea (2) se ejecutarán 3 OE (una asignación, una resta y una 
comparación) en cada una de las iteraciones del bucle más otras 3 al final, 
cuando se efectúa la salida del FOR. 

— Igual ocurre con la línea (3), también con 3 OE (una asignación, una suma y una 
comparación) por iteración, más otras 3 al final del bucle. 

— Y también en la línea (4), esta vez con 2 OE (asignación y comparación) más 
las 2 adicionales de terminación del bucle. 


— Por último, la línea (5) supone 2 OE (un incremento y una asignación). 


Con esto, el bucle interno se ejecutará ¡ veces, el medio (n—¿) veces, y el bucle 
exterior (n—1) veces, lo que conlleva un tiempo de ejecución de: 
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Tels $99 3(0[La=0)2))+>] De 
12 Else E 3+(4 )+2)+ Jere Els «$ aj+5)e3)]es- 


1+($ 0+Q0+0+10=0)+3))e3= 


3 2 e 
¡ ¿9 +15n* +13n ur O 
6 3 6 6 


b) Como el tiempo de ejecución es un polinomio de grado 3, la complejidad del 
algoritmo es de orden O(n”). 


Solución al Problema 1.6 (O) 


Para comprobar que O(f) < O(g) en cada caso y que esa inclusión es estricta, 
basta ver que lim,_...(An)g(n)) = 0, pues todas las funciones son continuas y por 


tanto los límites existen. Por consiguiente, 
O(1) < O(logn) < O(n) c O(nlogn) € O(1?) c O(n*) c O(n*) € 02”) c O(a!). 


Solución al Problema 1.7 (O) 


a) Por la definición de O, sabemos que fe O(g) si y sólo si existen c¡ > O y ny tales 
que fín) < c,g(n) para todo n >n,. 
Análogamente, por la definición de 2 tenemos que ge Q(f) si y sólo si existen cz 
> 0 y n, tales que g(n) > cof[n) para todo n >n>. Por consiguiente, 
>>) Si fe O(g) basta tomar c>=1/c, y nn, para ver que g(nje Q(f). 
<=) Recíprocamente, si ge Q(f) basta tomar c¡=1/c, y n¡=n, para que fe O(g). 
Obsérvese que esto es posible pues c, y cz son ambos estrictamente mayores que 
cero, y por tanto poseen inverso. 


po sin es par. > 
b) Sean An)> A E (Os 
l  sinesimpar. 

Entonces O(g)= om ), y por otro lado O(f)=0(1?), con lo cual f e O(n)=0(g). 
Sin embargo, si n es impar no puede existir c>0 tal que [(n) = 12 cn? = cg(n), y por 
consiguiente fe Q(g). 
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Intuitivamente, lo que buscamos es una función f cuyo crecimiento asintótico 
estuviera acotado superiormente por g (es decir, que f no creciera “más deprisa” 
que g) y que sin embargo fno estuviera acotado inferiormente por g. 

c) Veamos que log,ne O(log,n). 


Sabemos por las propiedades de los logaritmos que si a y b son números reales 
mayores que 1 se cumple que 


Pana log, n 


log, b ' 
Con esto, 


A a o 


n>0 log, n  n>w 


que es una constante real finita distinta de cero (pues a,b>1), y por tanto 
O(og,n) = O(log,n). 


Solución al Problema 1.8 (O) 


Para justificar estas afirmaciones nos apoyaremos en la definición de O y Q, y 
trataremos de encontrar la constante real c y el número natural ny que caracterizan 
las inecuaciones que definen a ambas cotas. 


— fi 0(3) pues n?< n?+ 1000n para todo n (podemos tomar c= 1, 19= 1). 
—fi£0() pues si n es impar no existen c y ny tales que n?< cn. 
— fis O(£) pues n*< n* si n> 100 (podemos tomar c= 1, ny= 101). 


— f.€ O(f1) pues basta tomar c y ny tales que c > 1 + 1000/n para todo n > ny que 
sabemos que existen pues (1000/n) tiende a cero. 


—fi£ 0(f) pues si n es impar no existen c y ny tales que n?+ 1000n < cn. 

- fe 0(£) pues n?+ 1000n < n* si n > 100 (podemos tomar c= 1, no= 101). 
—A£O(f) pues si n es par no existen c y ny tales que n*< cn?. 

—fi£O(f) pues si n es par no existen c y ny tales que n*< c(n?+ 1000n). 
fi 0() pues fin) < n”=fi(n) si n> 100 (podemos tomar c= 1, n9= 101). 
— fi£ O(A1) pues si n > 100 no existen c y no tales que n*< cn?. 

— fi O(f) pues si n > 100 no existen c y ny tales que n* < c(n?+ 1000n). 

— fu£ O(f) pues si n> 100, n impar, no existen c y no tales que n*< cn. 


fe Q(A) pues n> c(n?+ 1000n) para c y ny tales que c > 1 + 1000/n para todo 
n 2 ny que sabemos que existen pues (1000/n) tiende a cero. 


fit Q(f) pues si n es par no existen c y ny tales que n?> cn'. 
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— fit Q(f) pues si n> 100 no existen c y ny tales que n?> cn?. 
- fe Q(f) pues n”+ 100n > n? para todo n (podemos tomar c= 1, n= 1). 


— fEQ(% ues si n es par no existen c y no tales que n?+ 1000n > cn. 
p p y q 


—fe£Q(£) pues si n > 100 no existen c y ny tales que n?+ 1000n > cn”. 

— fi Q(f) pues si n es impar no existen c y ny tales que n > cn?. 

— fi Q(f) pues si n es impar no existen c y ny tales que n > c(n”+ 1000n). 
— fig Q(fi) pues si n es impar no existen c y ny tales que n > cn. 

— fis Q(A) pues n > mn? si n> 100 (podemos tomar c= 1, ny= 101). 

- fi Q(f) pues "> n?+ 1000n si n > 100 (podemos tomar c=1,n9)= 101). 
fi Q(fñ) pues > f(n) sin > 100 (podemos tomar c= 1, no= 101). 


Obsérvese que O(f,) = 01) = O(1?), O(4) = O(n?), Q(A) = Q(n), y O(5) = O(n). 


Solución al Problema 1.9 (9/8) 


Para resolver estas ecuaciones seguiremos generalmente el mismo método, basado 
en los resultados expuestos al comienzo del capítulo. Primero intentaremos 
transformar la ecuación en recurrencia en una ecuación de la forma: 


a, (n)+aT(n—-1)+..+a,T(n—-k)=bip,(n)+b)p,(n)+...+bp (n), 


para después resolver su ecuación característica asociada. Con las raíces de esta 
ecuación es fácil ya obtener el término general de la función buscada. 
a) T(n) = 3T(n-1) + 4T(n2) si n>1; 7(0)=0; 7(1) =1. 
Escribiendo la ecuación de otra forma: 
T(n) -3T(n-1) - 4T(n2) =0, 


ecuación en recurrencia homogénea con ecuación característica x—3x4 = 0, 
Resolviendo esta ecuación, sus raíces son 4 y —1, con lo cual: 


T(n) = c14" + ca 1). 


Para calcular las constantes necesitamos las condiciones iniciales: 


0=ST(D)=cF + (=D"=c +0 2 =1/5 
1=T()=c,4 +c,-1) =4e,-c,]  c,=-1/5 


Sustituyendo entonces en la ecuación anterior, obtenemos 
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Tm) == (4 (1 )eoW), 


b) T(n) = 27(n-1) — (n+5)3” si n>0; T(0) =0. 
Reescribiendo la ecuación obtenemos: 
T(n) -2T(n-1) = 4n+5)3”, 


que es una ecuación en recurrencia no homogénea cuya ecuación característica 
asociada es (x-2)(x-3)= 0, de raíces 2 y 3 (esta última con grado de multiplicidad 
dos), con lo cual 


T(n) = c12" + c23" + c3n3”. 


Para calcular las constantes necesitamos las condiciones iniciales. Como sólo 
disponemos de una y tenemos tres incógnitas, usamos la ecuación en recurrencia 
para obtener las otras dos: 


T(1) =27(0) -6:3=-18 
TQy=21(1)=7-9==99 


Con esto: 


0=T(0)=c,2+c,3 +c,0-3 =c,+c, c,=9 
-18=7()=c,2'+c,3 +c,1-3' =2c,+3c,+3c, +>c,=-9 
-99=T(2)=c,2* +0,3" +c,2-3” =4c, +9c, +18c, c¿=-3 


Sustituyendo estos valores en la ecuación anterior, obtenemos 


T(n) = 9-2" — 9-3" 3n3" eO(n3”). 


Existe otra forma de resolver este tipo de problemas, que se basa en manipular 

la ecuación original hasta convertirla en homogénea. Partiendo de la ecuación 7(n) 

2T(n-1) = (n+5)3", necesitamos escribir un sistema de ecuaciones basado en ella 
que permita anular el término no dependiente de 7(1). Para ello: 


— primero escribimos la recurrencia original, 
— la segunda ecuación se obtiene reemplazando n por n—1 y multiplicando por -6, 
— y la tercera se obtiene reemplazando n por n-2 y multiplicando por 9. 


De esta forma obtenemos: 


T(n) -2T(n-1) =  (n+5y3" 
-6T(n-1) +127(n2) =-6(n+4)3" 
9T(n2) —18T(n-3) = An+3)3"? 


Sumando estas tres ecuaciones conseguimos una ecuación homogénea: 


T(m) - 8T(n-1) + 21T(n-2) - 187(n-3) = 0 
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cuya ecuación característica es x* — 8x? + 21x — 18 = (1-23). A partir de aquí 
se puede resolver mediante el proceso descrito anteriormente. 


Como puede observarse, este método es más intuitivo pero menos metódico y 
ordenado que el que hemos utilizado para solucionar ecuaciones en recurrencia. 
Además, no hay una única forma de plantear esta ecuaciones. 


c) T(n) = 4T(1/2) + n? si n>4, n potencia de 2; T(1) = 1; TQ) =8. 
Haciendo el cambio n =2* (o, lo que es igual, k =logn) obtenemos 
TP) =4T04)) + 2% 
Llamando £;= T (2%, la ecuación final es 
(SA A, 
ecuación no homogénea con ecuación característica (x-4)?= 0. Por tanto, 
ti= 014 + coka! 
Necesitamos ahora deshacer los cambios hechos. Primero t;= 7(2*), con lo que 
TQ* = 0,4 + cok4'=c12% + crk27 


y después n =2*, obteniendo finalmente 


T(n) = có + con"logn. 

Para calcular las constantes necesitamos las condiciones iniciales: 
1=T()=c,P +c,1-0=C, el 
8=7(Q2)=c,2?+c,2* -1=4c, E Cs S 

Sustituyendo estos valores en la ecuación anterior, obtenemos 

T(n) = 1? + n'logn e O(n?logn). 
Existe otra forma de resolver este tipo de problemas, mediante el desarrollo 


“telescópico” de la ecuación en recurrencia. Eseribiremos la ecuación hasta llegar a 
una expresión en donde sólo aparezcan las condiciones iniciales: 


2 
Tím=4T|Llen=4dar Ll 2] +? = 
2 4) 2 
n 2 
=41T|%l+2=4 41% + 2] [+2 = 
4 8 4 


ae eo o =4T(1)+xm?. 
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De esta forma hemos ido desarrollando los términos de esta sucesión, cada uno 
en función de términos anteriores. Sólo nos queda por calcular el número de 
términos (x) que hemos tenido que desarrollar. 


Pero ese número x coincide con el número de términos de la sucesión 1/2, n/4, 
n/8,...,4,2,1, que es logn pues n es una potencia de 2. En consecuencia, 


T(n)= 4% T(D)+logn-n* =n'** 1+1l0g n-n? =n” +1l0g n- n* e O(1'logn). 


d) T(n) = 27(n/2) + nlogn si n>1, n potencia de 2. 
Haciendo el cambio n =2* (o, lo que es igual, k =logn) obtenemos 
TB =2704) + 10%. 
Llamando t¿= 7(24), la ecuación final es 
£i= Zty1 +12", 


ecuación en recurrencia no homogénea con ecuación característica (x-2)= 0. Por 
tanto, 


ti=02 + 040 F+ e ié2t 
Necesitamos ahora deshacer los cambios hechos. Primero t;= 7(2%), con lo que 
TP = 012% + c+ ct2*, 
y después n =2* (k=1logn), por lo cual 
T(n) = cin + canlogn + canlog?n. 


De esta ecuación no conocemos condiciones iniciales para calcular todas las 
constantes, pero sí es posible intentar fijar alguna de ellas. Para eso, basta sustituir 
la expresión que hemos encontrado para 7(n) en la ecuación original: 


nlogn = T(n) - 2T(n/2) = (cz — c,)n + 2c3nlogn, 
por lo que c3=c, y 2c3=1, de donde 
T(n) = cin + 1/2nlogn + 1/2nlog?n. 


En consecuencia 7(1)e O(nlog?n) independientemente de las condiciones iniciales. 


e) T(n) = 37(n/2) + 5n + 3 si n>1, n potencia de 2. 
Haciendo el cambio n =2* (o, lo que es igual, k=logn) obtenemos 
125 =3704%) +5:2%+3. 


Llamando £;= T (2%, la ecuación final es: 
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t= Mi + 5-24 +3, 


ecuación en recurrencia no homogénea cuya ecuación característica asociada es 
(1-3)(-2)0-1) = 0. Por tanto, 


== c¡3* + a + C3. 
Necesitamos ahora deshacer los cambios hechos. Primero t;= 7(2%), con lo que 
TP) =c13+ 072! + e, 
y después n =2* (k= logn), por lo cual 


1 log3 
T(n) = c13%%"+cm+c3=c nn" + can + Ca. 


De esta ecuación no conocemos condiciones iniciales para calcular todas las 
constantes, pero sí es posible intentar fijar alguna de ellas. Para eso basta sustituir 
la expresión que hemos encontrado para T(n) en la ecuación en recurrencia 
original, y obtenemos: 


cin + cn + a= 3(c1(18/3) + con/2 +3) + 5n +3. 

Igualando los coeficientes de n'“£ 
c3=-3/2 y c,—10, de donde 

T(n) = cin*É — 10n — 3/2. 


log3 


, n y los términos independientes obtenemos 


Como log3 > 1, 7(n) será de complejidad O(n 
T(n) e O(n) si c¡=0. 

Para ver cuándo c; vale cero estudiaremos los valores de las condiciones 
iniciales que le hacen tomar ese valor, en este caso 7(1). Por un lado, utilizando la 
ecuación original, tenemos que para n= 2: 


TQ)=3T(1) + 10+3. 
Por otro lado, basándonos en la ecuación que hemos obtenido, 
TQ) = 3c, -20 - 3/2. 


Igualando ambas ecuaciones, obtenemos que c; = T(1) + 23/2. Por tanto, 


) si c, es distinto de cero, o bien 


On) siT(1)+-23/2 
raya (207) TO) 
O(n) si 7 (1) = -23/2 
f) Tn) = 2T(1/2) + logn si n>1, n potencia de 2. 
Haciendo el cambio n =2* (o, lo que es igual, k =logn) obtenemos 
7125 =270*) + k. 


Llamando t¿= 7(24), la ecuación final es 
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£:= 2fj 1 + K, 
ecuación en recurrencia no homogénea que puede ser expresada como 

ti— 2t11=Kk 
y cuya ecuación característica asociada es (x-2Xx-1)?= 0. Por tanto, 

tp=012%+ c, + cak. 
Necesitamos ahora deshacer los cambios hechos. Primero 1¿= 7(2/), con lo que 
T2%) =4124+ 0, + 0k 
y después n =2* (k=10gn), y por tanto 
T(n) = cin + c, + calogn. 


De esta ecuación no conocemos condiciones iniciales para calcular todas las 
constantes, pero sí es posible intentar fijar alguna de ellas. Para eso, basta sustituir 
la expresión que hemos encontrado para T(n) en la ecuación en recurrencia 
original, y obtenemos: 


cin + c, + calogn =2(c,n/2 + c, + calogn — c3) + logn. 


Igualando los coeficientes de logn y los términos independientes obtenemos que 
c=—1 y c,=2, de donde 


T(n) = c¡n — 2 —logn. 


Esta función será de orden de complejidad O(n) si c, es distinto de cero, o bien 
T(n) e O(logn) si c,=0. 

Para ver cuándo c; vale cero estudiaremos los valores de las condiciones 
iniciales que le hacen tomar ese valor, en este caso 7(1). Por un lado, utilizando la 
ecuación original, tenemos que para n= 2: 


TQ) =27(1)+ 1. 
Por otro lado, basándonos en la ecuación que hemos obtenido 
TQ) =2c, -2-1. 
Igualando ambas ecuaciones, obtenemos que c; = T(1) +2. Por tanto, 
O(logn) siT(1) =-2 
T(me (logn) (1) 
O(n) s17(1) 4-2 
Es 1/2 02*. E 
g) T(n) =2T(n"") + logn con n=2* ; TQ)1. 
Haciendo el cambio n = 22 (k = loglogn) obtenemos la ecuación 


To02)=> 122 yuloga2” 
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k 
Llamando t¿= T (22 ), la ecuación final es 
t4= Ztia1 +2", 


ecuación en recurrencia no homogénea cuya ecuación característica es (x-2) =0. 
Por tanto, 


tl c12* + exar. 
Necesitamos ahora deshacer los cambios hechos. Primero t¿= TÍ 22%, con lo que 
k 
TQ2) =c/2% + cak2* 
k 
y después n = 22 (k= loglogn, o bien logn = 2*), por lo cual tenemos que 
T(n) = cilogn + calogn:loglogn. 


Para calcular las constantes necesitamos las condiciones iniciales. Como 
disponemos de sólo una y tenemos dos incógnitas, usamos la ecuación original 
para obtener la otra: 


T(4) = 272) + log4 = 4. 
Con esto: 
1=7(2)=c, log2+c, log2-0=c, c,=1 
4=T(4)=c,log4+c, log 4 -loglog4 =2c, +2c, 
Sustituyendo estos valores en la ecuación anterior, obtenemos 


T(n) = logn + logn:loglogn e O(logn:loglogn). 


h) T(n) = 57(1/2) + (nlognY si n>1, n potencia de 2; 7(1)=1. 
Haciendo el cambio n = 2* (o, lo que es igual, k=logn) obtenemos 
TQ =5TQ) + (0 1Y=5TQ%) + 1041. 
Llamando t¿= 7(24), la ecuación final es 
= ty + CAS, 
ecuación en recurrencia no homogénea que puede ser expresada como 
4-5 ti1= 4, 
cuya ecuación característica asociada es (x-5)(x-4)'= 0. Por tanto, 


=a5t J cr4* J cska* J calés!. 


Necesitamos ahora deshacer los cambios hechos. Primero t;= 7(2%), con lo que 
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Tr) = c¡5* c4* cska* cale4* 
y después n = ye (k = logn), por tanto 
T(m) = 5 + c,4" + cslogn4'*" + calog?n4 Y" = 


log5 log4 1 
cane + en" + cslognn 


og4 log4 — 


+ calog?n:n 


log5 2 2 2 2 
cin + con” + canlogn + canlogn. 


Para calcular las constantes necesitamos las condiciones iniciales. Como sólo 
disponemos de una y tenemos tres incógnitas, usamos la ecuación en recurrencia 
para obtener las otras dos: 

T(2) = 5T(1) + 2?=9; 
T(4) = 572) + 8?= 109; 
T(8) = 57(4) + 24?=1121. 


Con esto: 

1=7(1)=c, 1+c,1+c,1-0+c,1-0=c,+c, 
9=T(2)=c,5+c,4+c,4-1+c,4:1=5c, +4c, +4c, +4c, 

109 =T(4)=c,25+c,16+c,16-2+c,16-4=25c, +16c, + 32c, + 64c, 
1121=7(8) =c,125+c,64+c,64-3+c,64:-9=125c, + 64c, +192c, + 576c, 


Solucionando el sistema de ecuaciones obtenemos los valores de las constantes: 


c, =181 
c, = -180 
Cc, = 40 
c¿=-4 


y sustituyéndolos en la ecuación anterior, obtenemos 


T(n) = 18118 — 180n? — 40n"logn — 4n*log?n e Oprélog7n). 


i) Tm) = T(n-1) +27(n-2) -2T(n-3) si n>2; T(n) = 9n? — 15n + 106 si n=0,1,2. 


Reescribiendo la ecuación: 


T(n) - Tn-1) -2T(n2) + 27(n-3) =0, 
ecuación en recurrencia homogénea cuya ecuación característica asociada es 


—r-21+2=0. 


LA COMPLEJIDAD DE LOS ALGORITMOS 41 


Resolviendo esa ecuación, sus raíces son 1, y2 y—vY2, con lo cual 


Tm) = 0 +exd2 Y + es W2 y. 


Para calcular las constantes necesitamos las condiciones iniciales: 


106=T(0)=c, +c, +C, c, =100 
100=7()=c,+c,42-c,421>c,=3 ¿=> T(m)=100+ 312" (1+(-1)) 
112=7T(2)=c,+2c,+2c, c,=3 


En consecuencia, 7(n) e OQ”?). 
Y Tm =642)70/21/2T(1/41/n) si n>2, n potencia de 2;7(1)=1; 72)=3/2. 
Haciendo el cambio n = 2* (o, lo que es igual, k= logn) obtenemos 
1(25 = 8/DT2%) - (12724) - (1/2). 
Llamando t¿= 7(24), la ecuación final es 
ti= (3/2Mi1 — (1/2)Mro - (1/2) 

ecuación en recurrencia no homogénea en la forma 

(3/2) + (1/D6=-1/2Y, 
cuya ecuación característica asociada es (x-1)(x-1/2)= 0. Por tanto, 

£=c +02 ++ 
Necesitamos ahora deshacer los cambios hechos. Primero t;= 7(2%), con lo que 
TQ*) =01 +02 +02 * 

después n =2* (k =logn), con lo cual 

T(n) =c, + (1/n)c, + cx(logn/n). 


Para calcular las constantes necesitamos las condiciones iniciales. Como sólo 
disponemos de dos y tenemos tres incógnitas, usamos la ecuación en recurrencia 
para obtener la tercera: 


T(4) = B/D70) - U2AITO) - (1/4) = 3/2. 
Con esto: 
l1=T(1)=c, +c, c,=1 
312=T1(2)=06 +0 12) +c (12) +=, =0 A E ea e0(). 
3/2=T(4)=c, +c,(1/4)+c,(1/2) c,=1 Él 
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k) T(n) = 2T(1/4) + n'? si n>4, n potencia de 4. 
Haciendo el cambio n =2* (o, lo que es igual, k =logn) obtenemos 
1521100 
Llamando £;= T (2%, la ecuación final es 
t= 2£y 9 +2, 
ecuación en recurrencia no homogénea de la forma 


ti 2tp2= EAT 


EOS da ; 2 . 
cuya ecuación característica asociada es (x-2)(x-V2 ) = 0, o lo que es igual, 


(x+ 42 e 42 Y =0. Por tanto, 
a d2 Y + cd 2) + cs 2 Y. 


Necesitamos ahora deshacer los cambios hechos. Primero t¿= TÍ (29, con lo que 


T2%) = co V2 + e 2 + e 2 y 
después n = ar (k = logn), y por tanto obtenemos: 
Tn) =/n (a("E + ca + calogn). 


Si n es múltiplo de 4 entonces logn es par, y por tanto 1) 
Esto nos permite afirmar, llamando cy = c, + c,, que 


1" vale siempre 1. 
T(m)=4/n (co + calogn). 


De esta ecuación no conocemos condiciones iniciales para calcular todas las 
constantes, pero sí es posible intentar fijar alguna de ellas. Para eso, basta sustituir 
la expresión que hemos encontrado para 7(n) en la ecuación original: 


Jn (co + c3logn) = 2( Án 12 (co + calogn — 2c3)) + Yn. 


Igualando los coeficientes de Jn , Jn logn y los términos independientes 
obtenemos cz=1/2, de donde 
T(n) =/n (cy + 1/2logn) € e(V/n logn). 


Este problema también podría haberse solucionado mediante otro cambio, 
n= 4" (o, lo que es igual, k= log4n) obteniendo la ecuación 


T(4) =27(4*) + 4”. 


Esta nos lleva, tras llamar 1;= T (44, a la ecuación final 
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l= 2x1 + 2% 


cuya resolución conduce a la misma solución que la obtenida mediante el primer 
cambio. 


1) T(n) = 47(1/3) + n? si n>3, n potencia de 3. 
Haciendo el cambio n = 3* (o, lo que es igual, k =1log3n) obtenemos que 
T(3% =418%) +9, 
Llamando £;= T (32, la ecuación final es 
A 
ecuación en recurrencia no homogénea de la forma 
4 4t1=9, 
cuya ecuación característica asociada es (x-4)(x-9) = 0, Por tanto, 
ti= ci4* + cr. 
Necesitamos ahora deshacer los cambios hechos. Primero t;= 7(3%), con lo que 
T(3%) =014 + e,3% 
y después n = 3% (k= log3n), y por tanto 


1 2 logz 4 2 
Tra)=c PER =an e ECHA, 


De esta ecuación no conocemos condiciones iniciales para calcular todas las 
constantes, pero sí es posible intentar fijar alguna de ellas. Para eso, basta sustituir 
la expresión que hemos encontrado para T(n) en la ecuación en recurrencia 
original, y obtenemos: 


4 logz 4 % 2 4 
CS E a(3 +03) +n? cet Zo, +) 


Est y de n? obtenemos c,= 9/5, de donde 


Igualando los coeficientes de 1.” 
9 
T(n)=c,n** + o 


Como log,4 < 2, entonces T(n)e O(n?). 


Solución al Problema 1.10 (8) 
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a) Cierto. Se deduce de la propiedad 6 del apartado 1.3.1, pero veamos una posible 
demostración directa: 

Si 7,8 O(f), sabemos que existen c; > 0 y n; tales que T¡(n) € c¡f(n) para todo 
n > 1n;. Análogamente, como 7, O(f), existen cz > 0 y n, tales que Tr(n) < cam) 
para n 21». [1.1] 

Para comprobar que 7, + 7, € O(f), debemos encontrar una constante real c > 0 y 
un número natural ny tales que 7, (1) + Tan) < cf(n) para todo n > no. 

Apoyándonos en [1.1], basta tomar ny = máx(n,,n2y y c= c+ c,, con las que se 
verifica la ecuación [1.2] para todo n > no. 

Existe otra forma de demostrarlo, utilizando límites en caso de que estos 
existan, como sucede por ejemplo cuando las funciones son continuas: 


T,(m) 


Si T¡€ O(f), entonces lim ——-= k¡<oo, 
n= f(n) 
T 
Análogamente, como Te O(f), lim IS 01). ko, <oo, [1.3] 
as 10) 
T(M+T 
Veamos entonces que lim A. k<oo, [1.4] 


es — E) 


Pero [1.4] es cierto pues, como los dos límites en [1.3] son finitos y positivos 
podemos conmutar la suma con el límite y obtenemos que 


pim TADA) _ yy, TO py, TOD) 


a 1) A) FA) 


= k¡ + k,<oo, 


b) Cierto. 


Análogamente a lo realizado en el apartado anterior, si 7,€0(f), entonces 
. Tin Dn 
O (0) 


=k¡ < co. Igualmente, como Te O(f), lim k,< oo, [1.5] 
> f(n) n= f(n) 
T(n)-T 
Veamos entonces que lim EE). k<oo, [1.6] 


> 
Pero [1.6] es cierto pues, como los dos límites en [1.5] existen y son finitos y 
positivos podemos conmutar la resta con el límite y obtenemos que 
T.(n)=T EE TL 
im LG) TD) _ y, LD) — y, DD 
E ¿FM sl. 8 500) 


=k¡-k,<oo, 
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c) Falso. 


Consideremos T;(n) = n?, Tn) =1n, y fin) = n?. Tenemos por tanto que 7, O(f) 
y Te O(f), pero sin embargo T ¡(1 /T(n) =n2£ 0(1). 


d) Falso. 


Consideremos de nuevo T;¡(n) = 17, TA(n) =1, y Án) = n?. Tenemos por tanto que 
T¡e O(f) y T2e O(f), pero sin embargo 7,2 O(T,) pues ng O(n). 


Solución al Problema 1.11 (O) 


n?, sines par. 
Sean fl) =n y g(n) = AS 
l,  sinesimpar. 
S1 n es impar, no podemos encontrar ninguna constante c tal que 
An) =n<cg(n) =cC, 


y por tanto f20O(g). Por otro lado, si n es par no podemos encontrar ninguna 
constante c tal que 


g(n)=1"<cfn) =cn, 
y por tanto gé O(f). 


Solución al Problema 1.12 (O) 


Para comprobar que log'ne O(n) basta ver que 


 logn 
lim 8% =q 

n—>0 n 
para todo k. Pero eso es cierto siempre. Obsérvese además que por esa misma razón 
log'ng Q(n) para cualquier k> 0, 


Solución al Problema 1.13 


Vamos a suponer que los tiempos de ejecución de las funciones Esvacio, [zq, Der y 
Raiz es de c operaciones elementales (OE), que el tiempo de ejecución de Opera es 
d OE, y el de Max2 es 1 OE. 


Procedimiento /Inorden (O) 
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Para calcular el tiempo de ejecución, calcularemos primero el número de 
operaciones elementales (OE) que se realizan: 


En la línea (1) se ejecutan 2+c OE: la llamada a Esvacio (1 OE), el tiempo de 
ejecución de este procedimiento (c) y una negación. 

En la línea (2) se efectúa la llamada a /zq (1 OE), lo que tarda ésta en ejecutarse 
(c OE) más la llamada a [norden (1 OE) y lo que tarde ésta en ejecutarse, que va 
a depender del número de elementos del árbol /zq(1). 

En la línea (3) se ejecutan 2+c+d OE: dos llamadas a procedimientos y sus 
respectivos tiempos de ejecución. 


El número de OE de la línea (4) se calcula de forma análoga a la línea (2): 2+c 
más lo que tarda Inorden en ejecutarse con el número de elementos que hay en 
Der(t). 


Para estudiar el tiempo de ejecución, vamos a considerar dos casos extremos: 


que el árbol sea degenerado (es decir, una lista) y que sea equilibrado. Cualquier 
árbol se encuentra en una situación intermedia a estos dos casos. 


Si £ es degenerado, podemos suponer sin pérdida de generalidad que 
Esvacio(Izq(t)) y que para todo a subárbol de £ se verifica que Esvacio(Izq(a)). 
Por tanto, el número de OE que se realizan en la ejecución de Inorden(t) para un 
árbol £ con n elementos es: 


T(n) = (Q+c) + Q+c+T(0)) + Q+c+d) + Q+c+T(n-1)) = 
8+4c+d+T(0)+T(n-1). 
T(0) = 2+c. 


Con esto, T(n) = 10 + 5c + d + T(n—1), ecuación en recurrencia no homogénea 
que podemos resolver desarrollándola telescópicamente: 


T(n) =10+5c+ d+ T(n-1)= (10 +5c+d)+(10+5c+d)+ T(n2)=...= 
(10 + Sc + din + Q+c) e O(n) 


Si £ es equilibrado sus dos subárboles (izquierdo y derecho) tienen del orden de 
n/2 elementos y son a su vez equilibrados. Por tanto, el número de OE que se 
realizan en la ejecución de Inorden(t) para un árbol £ con n elementos es: 


T(nyF Q+c) + Q+c+T(n/2)) + Q+c+d) + Q+c+T(n/2)) = 8+4c+d+2T(n/2). 
T(0F 2+c. 


eN . . k 
Para resolver esta ecuación en recurrencia se hace el cambio 1¿= 7(25), con lo 
que obtenemos 


UNA 2tr1=8 E Ac 4 d, 


ecuación no homogénea con ecuación característica (x-2)(x-1) = 0. Por tanto, 
e k 
ti=C12+C) 
y, deshaciendo los cambios, 


T(n) =cin +C>. 
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Para calcular las constantes, nos apoyamos en la condición inicial 7(0)=2+c, 
junto con el valor de 7(1), que puede ser calculado basándonos en la expresión 
de la ecuación en recurrencia: 7(1) =8 + 4c + d+ 2(Q2 + c), obteniendo 


T(n) = (10 + 5c+ din + (2+c) e O(n). 


Función Altura (9) 


Para determinar el tiempo de ejecución, calcularemos primero el número de 
operaciones elementales (OE) que se realizan: 


— En la línea (1) se ejecutan 1+c OE: la llamada a Esvacio (1 OE) más el tiempo 
de ejecución de este procedimiento (c OE). 


— En la línea (2) se realiza 1 OE. 
— Enla línea (4) se efectúan: 


a) la llamada a /zq (1 OB), lo que tarda ésta en ejecutarse (c OE) más la llamada 
a Altura (1 OE) y lo que tarde ésta en ejecutarse, que va a depender del 
número de elementos del árbol /zq(1); más 


b) la llamada a Der (1 OE), lo que tarda ésta en ejecutarse (c OE) más la 
llamada a Altura (1 OE) y lo que tarde ésta en ejecutarse, que va a depender 
del número de elementos del árbol Der(t); más 


c) el cálculo del máximo de ambos números (1 OE), un incremento (1 OE) y el 
RETURN (1 OB). 


Para estudiar el tiempo de ejecución de esta función consideraremos los mismos 
casos que para la función /norden: que el árbol sea degenerado (es decir, una lista) 
o que sea equilibrado. 


e Si f es degenerado, podemos suponer sin pérdida de generalidad que 
Esvacio(lzq(t)) y que para todo a subárbol de £ se verifica que Esvacio(Izq(a)). 
Por tanto, el número de OE que se realizan en la ejecución de Altura(t) para un 
árbol t con n elementos es: 


T(M= (1+c) + (1+c+1+T(0) + 1+c+1+T7(n-1)) +3 =8 + 3c + T(0) + T(n-1). 
T(0 (1+c) + 1 =2+c. 


Con esto, 7(n)=10+4c+T(n-1), ecuación en recurrencia no homogénea que 
podemos resolver desarrollándola telescópicamente: 


T(n) = 10 + 4c + T(n-1) = (10 + 4c) + (10 + 4c) + T(n2)=...= 
(10 +40n+(Q+c)€ O(n) 


e Sites equilibrado sus dos subárboles tienen del orden de n/2 elementos y son 
también equilibrados. Por tanto, el número de OE que se realizan en la 
ejecución de A/tura(t) para un árbol £ con n elementos es: 


T(nyF (1+c) + (1+c+1+7(n/2)) + 1+c+1+T(1/2)) +3 = 8 + 3c + 27(n/2). 
T(0F 2+c. 
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¿7 . : k 
Para resolver esta ecuación en recurrencia se hace el cambio 1¿=7(25), con lo que 
obtenemos 


ti— 2f11= 8 + 3c, 


ecuación no homogénea de ecuación característica (x-2)(—1) = 0. Por tanto, 


t1= 02 C2. 


Deshaciendo los cambios, 


T(n) = cin + Cc). 


Para calcular las constantes, nos apoyamos en la condición inicial 


T(0)2+c, junto con el valor de 7(1), que puede ser calculado basándonos en la 
expresión de la ecuación en recurrencia: 7(1) = 8 + 3c + 2(2 + c). Finalmente 
obtenemos 


T(n) =(10+40n +(Q+c) € O(n). 


Función Mezcla (8) 


Para resolver este problema vamos a suponer que el tiempo de ejecución del 
procedimiento /ns, que inserta un elemento en un árbol binario de búsqueda, es 
Alogn+B, siendo A y B dos constantes. Supongamos también que n y m son el 
número de elementos de £/ y £2 respectivamente. 


Para estudiar el tiempo de ejecución 7(m,m) consideraremos, al igual que 
hicimos para la función anterior, dos casos extremos: que el árbol 12 sea 
degenerado (es decir, una lista) o que sea equilibrado. 


+ Si t2 es degenerado, podemos suponer sin pérdida de generalidad que 
Esvacio(Izq(t2)) y que para todo a subárbol de 12 se verifica que 
Esvacio(Izq(a)). Por tanto, vamos a ver el número de OE que se realizan en cada 
línea de la función en este caso: 


En la línea (1) se invoca a Esvacio(t1), lo que supone 1+c OE. 
En la línea (2) se efectúa 1 OE. 
Análogamente, las líneas (3) y (4) realizan (1+c) y 1 respectivamente. 


Para estudiar el número de OE que realiza la línea (6), vamos a dividirla en 
cuatro partes: 


a) al:=Ins(t1,Raiz(t2)), siendo al una variable auxiliar para efectuar los 
cálculos. Se efectúan 2+c+4logn+B operaciones elementales: la llamada a 
Raiz (1), el tiempo que ésta tarda (c), la llamada a Ins (1 OE), y su tiempo 
de ejecución (Alogn+B). 

b) a2:=Mezcla(al,Izq(t2)), siendo a2 una variable auxiliar para efectuar los 
cálculos. Se efectúan aquí 2+c+T(n+1,0) operaciones elementales: 
llamada a Izg (1), el tiempo que ésta tarda (c), la llamada a Mezcla 
(1 OE), y su tiempo de ejecución, que será 7(n+1,0), pues estamos 
suponiendo que Esvacio(Izg(a)) para todo a subárbol de £2. 


c) a3:=Mezcla(a2,Der(t2)), siendo a3 una variable auxiliar para efectuar los 
cálculos. Se efectúan 2+c+T(nm+1,m-1) operaciones elementales: la 
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llamada a Der (1), el tiempo que ésta tarda (c), la llamada a Mezcla 
(1 OE), y su tiempo de ejecución, que será 7(n+1,m-1), pues estamos 
suponiendo que Esvacio(lzq(a)) para todo a subárbol de 12 o, lo que es 
igual, que el número de elementos de Der(t) es m-1. 


d) RETURN a3, que realiza 1 OE. 
Por tanto, la ejecución de Mezcla(t1,t2) en este caso es : 


T(n,m) =9 + 5c + B + Alogn + T(n+1,0) + T(n+1,m-1) 


con las condiciones iniciales 7(0,m) =2 +c y T(n,0) = 3 + 2c. Para resolver la 
ecuación en recurrencia podemos expresarla como: 


T(n,m) = 12 +7c + B + Alogn + T(n+1,m-1) 


haciendo uso de la segunda condición inicial. Desarrollando telescópicamente la 
ecuación: 


T(n,m) = 12 +7c + B + Alogn + T(n+1,m-1) = 
= (12+7c+B+4logn) + (12+7c+B+Alog(n+1)) + T(n+2,m-2) = 


m-1 


m(12+7c+B)+ (Sa log(n + 0) +T(n+m,0)- 


i=0 


m-1 
m(12+7c+ B)+2c+3+ e Host + 0) 
i=0 


Pero como log(n+i) <€ log(n+m) para todo 0<i¡<m, 


T(n,m) € m(2 + Tc + B) + 2c + 3 + Amlog(n+m) e O(mlog(n+m)) 


e El segundo caso es que £2 sea equilibrado, para el que se demuestra de forma 
análoga que 


T(n,m) e O(mlog(n+m)). 
Solución al Problema 1.14 (O) 


Para comprobar que O(f)CO(g), basta ver que lim,,.. (An)/g(n)) = 0 en cada caso 


pues las funciones son continuas, lo que implica la existencia de los límites. De 
esta forma se obtiene la siguiente ordenación: 


0((1/3)") < 0(17) < O(loglogn) < O(logn) < O(log?n) € O(/n ) € On log?n) 
c O(n/logn) c O(n) c O(n?) c 0((3/2)”. 


Solución al Problema 1.15 (8) 
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Para resolver la ecuación 


1 n-1 
T(n) = ($ 10) +cn, 
AX i=0 
siendo 7(0) = 0, podemos reescribirla como: 


n-1 
nT(n) = Y T(1)+cn? [1.7] 
i=0 


Por otro lado, para n—1 obtenemos: 


(n-—DT(n-1)= STi) +c(n—-1y [1.8] 


Restando [1.7] y [1.8]: 
nT(n) - nT(n-1) + T(n-1) = T(n-1) + cQn-1) > nT(n) = nT(n-1) + c(Qn-1) > 


T(n) = T(n-1) + cQ-1/n). 
Desarrollando telescópicamente la ecuación en recurrencia: 
T(n) = Tn-1)+cQ-1/n)= 


= T(n2) + cQ — 1/(n-1)) + cQ — 1/n) = 
= T(n3) + cQ — 1(m2)) + cQ — 1/(n-1)) + cQ — 1/n)= 


= r0o)+0d[2-;)- 
1 


i=1 


ya que teníamos que 7(0) = 0. Veamos cual es el orden de 7(m): 


a) Como (-1/i) < 2 para todo ¡> 0, T(n) < cd 2 = 2cn > T(nje O(m). 


i=l 
b) Como (2-1/i) > 1 para todo ¡>0, T(n)> cd 1 =cn => T(meQ(n). 


¡=l 
Por tanto, T(n)e O(n). 


Solución al Problema 1.16 
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Función BuscBin (O) 


a) Para determinar su tiempo de ejecución, calcularemos primero el número de 
operaciones elementales (OE) que se realizan: 


— En la línea (1) se ejecutan la comparación del /F (1 OE), y un acceso a un 
vector (1 OE), una comparación (1 OE) y un RETURN (1 OE) si la 
condición es verdadera. 


— En la línea (3) se realizan 3 OE (suma, división y asignación). 

— En la línea (4) hay un acceso a un vector (1 OE) y una comparación (1 OE), 
y además 1 OE en caso de que la condición del IF sea verdadera. 

— En la línea (5) hay un acceso a un vector (1 OE) y una comparación (1 OE). 


— Las líneas (6) y (8) efectúan 3+7(n/2) cada una: una operación aritmética 
(incremento o decremento de 1), una llamada a la función BuscBin (lo que 
supone 1 OE), más lo que tarde en ejecutarse la función con la mitad de los 
elementos y un RETURN (1 OB). 

Por tanto obtenemos la ecuación en recurrencia 7(n) = 11 + 7(n/2), con la 

condición inicial 7(1) = 4. Para resolverla, haciendo el cambio 1; = 7(2%) 

obtenemos 


ti tu=11, 
ecuación no homogénea cuya ecuación característica es (x-1)?= 0. Por tanto, 
t4= Cik + Cc) 
y, deshaciendo los cambios, 
T(n) = c¡logn + c>. 


Para calcular las constantes, nos basaremos en la condición inicial 7(1) = 4, 
junto con el valor de 7(2), que podemos calcular apoyándonos en la expresión 
de la ecuación en recurrencia: 7(2) = 11 + 4= 15. Finalmente obtenemos 


T(n) = 11logn + 4 e O(logn) 


b) La recursión de este programa, por tratarse de un caso de recursión de cola, 
puede ser eliminada mediante un bucle que simule las llamadas recursivas a la 
función. La condición de terminación del bucle puede ser tomada del caso base 
de la función recursiva y el cuerpo de dicho bucle consiste en una preparación 
de los argumentos de la función recursiva y el cálculo que ésta realiza: 


PROCEDURE BuscBlIt(a:vector;prim,ult:CARDINAL;x: INTEGER) : BOOLEAN; 
VAR mitad: CARDINAL; 


BEGIN 
WHILE (prim<ult) DO (* 1 1) 
mitad:=(prim+ult)DIV 2; (x* 2 x) 
IF x=a[mitad] THEN RETURN TRUE (* 3 *) 
ELSIF (x<a[mitad]) THEN (* 4 *) 
ult:=mitad-1 (x* 5 x) 
ELSE (* 6 *) 
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c) 
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prim:=mitad+1 (x* 7 *) 

END (x* 8 *) 
END; (x* 9 *) 
RETURN x=a [ult] (x* 10 x) 


END BuscBlt; 


Para el cálculo del tiempo de ejecución y la complejidad de la función no 
recursiva podemos seguir un proceso análogo al que seguimos para la función 
Algoritmo2 (en el problema 1.5). Para determinar el tiempo de ejecución, 
calcularemos primero el número de operaciones elementales (OE) que se 
realizan en cada una de las líneas: 


— En la línea (1) se efectúa la condición del bucle, que supone 1 OE (la 
comparación). 

— Las líneas (2) a (9) componen el cuerpo del bucle, y contabilizan 3, 2+1, 2, 
2,0, 2,0 y 0 OE respectivamente. 


— Por último, la línea (10) supone 3 OE. A ella se llega cuando la condición 
del bucle deja de verificarse. 


El bucle se repite hasta que su condición sea falsa, acabando la función al 
ejecutarse la línea (10). Cada iteración del bucle está compuesta por las líneas 
(1) a (9), junto con una ejecución adicional de la línea (1) que es la que ocasiona 
la salida del bucle. En cada iteración se reduce a la mitad los elementos a 
considerar, por lo que el bucle se repite logn veces. Por tanto, en el peor caso, 


logn 
T(n) = [(Ea+a+2+2+2)+1)+3- 10logn + 4e O(logn). 


i=l 


Como puede verse, el tiempo de ejecución de ambas funciones es 
prácticamente igual, lo que a priori implica que cualquiera de las dos pueden 
usarse indistintamente. Sin embargo, hay que tener en cuenta la mayor 
complejidad espacial que siempre suponen los procedimientos recursivos por la 
utilización de la pila, lo que hace que ante una igualdad de tiempos de 
ejecución, los procedimientos iterativos sean preferibles frente a los recursivos. 
Pero no sólo la complejidad ha de ser tenida en cuenta para la elección del 
algoritmo. La claridad y sencillez del código es un factor también a considerar, 
pues ello va a implicar una mejor legibilidad y una depuración del programa y 
mantenimiento más fácil, aspectos todos ellos muy importantes. 


Función Sumadigitos (8) 


a) 


Para calcular el tiempo de ejecución, calcularemos primero el número de 

operaciones elementales (OE) que se realizan: 

— En la línea (1) se ejecutan una comparación (1 OE) y un RETURN (1 OE) si 
la condición es verdadera. 

— En la línea (2) se efectúa una división (1 OE), una llamada a la función 


Sumadigitos (1 OE), más lo que tarda ésta con un décimo del tamaño de su 
entrada, una suma (1 OE), un resto (1 OE), y un RETURN (1 OE). 
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b) 


c) 


Llamando n al parámetro num de la función, obtenemos la ecuación en 
recurrencia 7(n) = 6 + 7(n/10), con la condición inicial 7(1) =2. 


Para resolverla hacemos los cambios n = 10% (o, lo que es igual, k= lo0g¡0n) y 
t¿= T(10*) y obtenemos 
li t4156, 
ecuación no homogénea cuya ecuación característica es (x-1)?= 0. Por tanto, 
t= C¡ + Cok. 
Deshaciendo los cambios, 
T(n) =c¡ + calog¡0n. 


Para calcular las constantes, nos apoyamos en la condición inicial 7(1) = 2, 
junto con el valor de 7(10), que puede ser calculado apoyándonos en la 
expresión de la ecuación en recurrencia: 7(10) = 6 + 2 = 8, Finalmente 
obtenemos 


T(n) = 6 log ion + 2 e O(logn) 


Como vemos, en esta caso la complejidad de la función depende del 
logaritmo en base 10 de su parámetro num (esto es, de su número de dígitos). 


La recursión de este algoritmo puede ser eliminada mediante un bucle que 
simule las llamadas recursivas a la función, cuya condición de terminación 
puede ser tomada del caso base de la función recursiva, y cuyo cuerpo consiste 
en los cálculos que ésta realiza, junto con una preparación de los argumentos de 
la siguiente llamada. En concreto, el algoritmo que implementa el algoritmo no 
recursivo es el siguiente: 


PROCEDURE Sumadigitos_it (num: CARDINAL) : CARDINAL; 
VAR s:CARDINAL; 


BEGIN 
s:=num MOD 10; (x* 1 1) 
WHILE num>=10 DO (* 2 *) 
num:=num DIV 10; (x* 3 x) 
s:=s+ (num MOD 10) (x* 4 x) 
END; (* 5 x*) 
RETURN s (* 6 *) 


END Sumadigitos_it; 


Para determinar el tiempo de ejecución, calcularemos primero el número de 
operaciones elementales (OE) que se realizan: 
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— En la línea (1) se ejecutan 2 OE (un resto y una asignación). 
— En la línea (2) se efectúa la condición del bucle, que supone 1 OE. 


— Las líneas (3) y (4) componen el cuerpo del bucle, y contabilizan 2 y 3 OE 
respectivamente. 


— Por último, la línea (6) supone 1 OE. A ella se llega cuando la condición del 
bucle deja de verificarse. 


Con esto, se efectúa la línea (1), y después se repite el bucle hasta que su 
condición sea falsa, acabando la función al ejecutarse la línea (6). Cada 
iteración del bucle está compuesta por las líneas (2) a (4), junto con una 
ejecución adicional de la línea (2) que es la que ocasiona la salida del bucle. En 
cada iteración se diezman los elementos a considerar, por lo que el bucle se 
repite log¡pn veces. 


log¡o 1 


Por tanto, T(n) =2 + ( ) (1+24+ 2) + ) +1=6log,,n+4 e O(logn). 


i=l 


Llegados a este punto vemos que ocurre lo mismo que en el algoritmo 
anterior. Los tiempos de ejecución de las funciones recursiva e iterativa son 
similares. Es por tanto una decisión del usuario decantarse por el diseño 
generalmente más robusto ofrecido por los algoritmos recursivos, frente a la 
menor complejidad espacial que presentan los iterativos al no utilizar la pila de 
recursión. 


Solución al Problema 1.17 (O) 


a) Para determinar el tiempo de ejecución del algoritmo, calcularemos primero el 
número de operaciones elementales (OE) que se realizan: 


En la línea (1) se puede ejecutar o sólo la condición (si ésta es falsa) con un 
total de 1 OE, o bien la sentencia entera, con un total de 3 OE. 


Las líneas (2) y (3) están compuestas por operaciones aritméticas y 
asignaciones, y se realizan un total de 3 OE en cada una. 


La línea (4) tiene tres partes: un acceso a a|mitad|, que supone 1 OE; la llamada 
a Raro(a,prim,prim+terc), que supone 2 + 7(n/3) (1 de la suma de prim y terc, 1 
OE de la llamada a Raro, y luego el tiempo de ejecución de Raro para un tercio 
del número de elementos con los que fue invocada la función original); y la 
llamada a Raro(a,ult-terc,ult), que supone 2 + T(n/3) OE por la misma razón. El 
resultado de las tres partes ha de sumarse (lo que supone 2 OE) y luego hacer un 
RETURN (1 OE). En resumen, en la línea (4) se ejecutan un total de 1 + 2 + 
T(n/3) + 2 + T(n/3)+2+1=8 +27(n/3) OE. 


Con esto: 


1. 
2. 


Sin= 1 (caso base), se ejecuta sólo la línea (1) y por tanto 7(1) = 3. 


Si n=3%*, con k> 0, se ejecuta sólo la condición del /F (1) y luego el resto de las 
líneas (2) a (4), con lo cual 7(n) = 15 + 27(113). 
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Tenemos por tanto el tiempo de ejecución del algoritmo definido mediante una 
ecuación en recurrencia. Para resolverla, haciendo el cambio n = 3% queda 


1(3%) =21(3%) + 15 
y llamando 1, a 7(3*) obtenemos 
1= ti +15, 


ecuación en recurrencia no homogénea de ecuación característica (x-2(x-1) = 0 y 
consecuentemente 


(1 Ear J iS Ez F Ca. 
Cambiando entonces 1; por 7(3*) queda 
T(3*) = e12+ C2, 


y deshaciendo el cambio n = 3* (o, lo que es igual, k = log3n), obtenemos 
finalmente 


T(n) = c12108 + (,= e nl0832 + c), 


Para calcular las constantes necesitamos resolver un sistema de dos ecuaciones 
con dos incógnitas (c, y c,) basándonos en dos condiciones iniciales de la ecuación 
en recurrencia. Como de partida sólo disponemos de una (7(1) = 3), necesitamos 
obtener otra. Para ello, apoyándonos en la definición recursiva de T y para n = 3, 
obtenemos 


T3)=15 +27(0/3)=15+27(1)=21. 


Por tanto 


3=T(1)=c,+<c, c, =18 
> 
21=7(3) =2c, +c, c, =-15 


Sustituyendo estos valores en la ecuación, obtenemos finalmente 


T(n) = 18nl082 — 15, 


b) T(n)e O(n!o832), 
Para justificarlo, basándonos en las propiedades de O, basta ver que 
T 
lim 07 
n>0 y 083 


existe, es acotado y distinto de cero. Pero eso es cierto ya que 


TÉCNICAS DE DISEÑO DE ALGORITMOS 


1 log; 2 =1 
lim LD) - q Sl =18. 
n>0 y 3 n>0 n 


